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Chapter 1 


Introduction 


The study of the Lee-Yang model is important for a general understanding of two di¬ 
mensional integrable models. The main motivation behind that comes from AdS/CFT 
duality. In our quest for understanding realistic but very complicated models like the 
Super Yang-Mills gauge theories, there is an important conjecture called the AdS/CFT 
duality and it jl] states the equivalence of A/" = 4 Super Yang-Mills gauge theory 
with superstrings on AdS^ x . The correspondence is extremely interesting as it 
links the very difficult non-perturbative physics of gauge theory to (semi) classical 
string/supergravity theory. 

As such it allows to gain new insight into various gauge theoretical phenomena but at 
the same time makes it very difficult to test and prove. A real breakthrough in this 
respect is the discovery of integrability on both sides of the duality la ism El Ell?]. On 
the string theory side it means that the light-cone quantized worldsheet sigma model 
is an integrable quantum field theory, while on the gauge theory side it manifests itself 
in the appearance of spin chains. 

It becomes obvious that we need to have a deeper understanding of two dimensional 
integrable quantum field theoretical models as a starting step to understand the more 
complicated theories leading to realistic models. 

In this thesis I choose the Lee-Yang model and go through different approaches to 
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analyze the model using the form factor approach and the bootstrap program, the 
lattice model and the TBA equations from the lattice as different approaches that 
lead to a full picture about the model. 

The bootstrap program aims to classify and explicitly solve 1+1 dimensional inte- 
grable quantum held theories by constructing all of their Wightman functions ( for a 
recent review [9] and references [ID1EIEI[I31[I511891E]). In the hrst step, called the 
S-matrix bootstrap, the scattering matrix, connecting asymptotic in and out states, 
is determined from its properties such as factorizability, unitarity, crossing symmetry 
and Yang-Baxter equation (YBE) supplemented by the maximal analyticity assump¬ 
tion |21]. In the second step, called the form factor bootstrap, matrix elements of local 
operators between asymptotic states are computed using their analytical properties 
originating from the already computed S-matrix. Supposing maximal analycity leads 
to a set of solutions each of which corresponds to a local operator of the theory |8]. 
In the third step these bulk form factors are used to build up the correlation (Wight¬ 
man) functions via their spectral representations and describe the theory completely 
off mass shell. This program has been implemented for a wide range of theories as in 

[iHiiiHiiniiiniEniEilEaEa]. 

The analogous bootstrap program for 1+1 dimensional integrable boundary quan¬ 
tum held theories has been already developed. The hrst step is called the R-matrix 
bootstrap j25]: In boundary theories the asymptotic states are connected by the R- 
matrix, which, as a consequence of integrability factorizes and satishes the unitarity, 
boundary crossing unitarity and the boundary YBE (BYBE) requirements. These 
equations supplemented by the maximal analytical assumptions makes it possible to 
determine the rehection matrices and provide the complete information about the the¬ 
ory on mass shell. In the second step we are interested in the matrix elements of local 
operators localized both in the bulk and also at the boundary. Due to the absence 
of translational invariance the bulk operators' one point functions acquire nontrivial 
space dependence which can be calculated in the crossed channel using the knowledge 
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of the boundary state together with the bulk form factors |43| . 

For the matrix elements of local boundary operators axioms can be derived from their 
analytical properties originating from the already computed R-matrix j29j. Supposing 
maximal analytical leads to a set of solutions each of which corresponds to a local 
boundary operator of the theory and is uniquely related to a vector in the ultravi¬ 
olet Hilbert space. The explicit form of the boundary form factors via the spectral 
representation of the boundary correlation functions provides a partial description of 
the theory off the mass shell. A full description would include correlation functions 
of operators localized in bulk as well, but this complicated problem has not been 
addressed yet. 

Since any two dimensional defect theory can be mapped to a boundary theory j31] the 
development of a separate bootstrap program for their solution seems to be redundant. 
However, integrable defect theories are severely restricted and one can go much beyond 
the boundary bootstrap program explained above: We can determine the form factors 
of both types of operators, those localized in the bulk and also the ones localized 
on the defect. With the help of these form factors we are able to derive spectral 
representation for any correlation function and in principle fully solve the theory off 
the mass shell as we show in |33) . 

In developing a defect form factor program the first step is the T-matrix bootstrap. 
Interacting integrable defect theories are purely transmitting [3D1 EH E2] and topo¬ 
logical. As a consequence a momentum like quantity is conserved [35l l36] and the 
location of the defect can be changed without affecting the spectrum of the theory 
[SZlEH]. This fact, together with integrability lead to the factorization of scattering 
amplitudes into the product of pairwise scattering and individual transmission and 
enable one to determine the transmission factors from defect YBE (DYBE), unitarity 
and defect crossing unitarity |39l HH 00]. The second step is the defect form factor 
bootstrap: Once the transmission factors are known we can formulate the axioms that 
have to be satisfied by the matrix elements of local defect operators. We will analyze 
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both operators localized in the bulk and also on the defect. By hnding their solutions 
the spectral representation of any correlator can be determined and theory can be 
solved completely. 

In |53l El ESI ESI ET], Thermodynamic Bethe Ansatz (TBA) equations have been 
introduced as an important tool in the study of both massive and massless integrable 
quantum held theories. Extensive studies have been carried out on scaling energies of 
vacuum or ground states. However only relatively few excited states jSSl El El EDI ED] 
were possible by TBA analysis and these are primarily restricted to massive and 
diagonal scattering theories. So despite considerable successes, the application of 
TBA methods was limited. The primary obstacle is that there is no systematic and 
unihed derivation of excited state TBA equations m 

The Lee-Yang model was studied from the TBA approach. The periodic Lee-Yang 
was analyzed in for the groundstate, and in isa and [90] for the excited states. 
The equations were solved based on assumptions about the analytic structure of the 
model, and were also supported by numerical results from the TCSA. 

The TBA approach was also used to study the boundary Lee-Yang in groundstate 
|91| and also in excited states, [901 where as the defect groundstate was analyzed in 

0Z|. 

However, the lattice approach is far more reaching. It is a systematic approach that 
allows to obtain both massive and massless excited TBA equations by studying the 
continuum scaling limit of the associated integrable lattice models. The most impor¬ 
tant input from the lattice approach is an insight into the analytic structure of the 
excited state solutions of the TBA equations. Previously this structure had to be 
guessed. In contrast, in the lattice approach, the analytical structure can be probed 
by direct numerical calculations of hnite size transfer matrices. 

The lattice model is very general and was used to study several models like the tri- 
critical Ising model HHilHI, by considering the massive tricritical Ising model Ad (4,5) 
perturbed by the thermal operator 
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There have been many relevant studies of the A 4 lattice model and the more general A^ 
models from the lattice viewpoint. For the ^4 model, the off-critical TEA functional 
equation for periodic boundary conditions has been derived and solved |M1169] for the 
bulk properties and correlation lengths. The off-critical TEA functional equations for 
the AL models were derived by Kliimper and Pearce izniEiiizs]. Eut only the critical 
or “conformal TEA" equations were derived and solved in the critical scaling limit for 
the central charges and conformal weights. The very same off-critical TEA functional 
equations for Al models were subsequently derived iza in the presence of integrable 
boundaries showing that the TEA functional equations are universal in the sense that 
they are independent of the boundary conditions. 

In this thesis we turn our attention to the simplest example of a non-unitary minimal 
theory Ai{p,p') with \p — p'\ 7 ^ 1, namely, the Lee-Yang minimal model Ad(2,5) j74| . 

Here we study the Lee-Yang model on the lattice. We analyze the periodic, boundary 
and the seam cases in both massive and massless regimes. We derive their ground 
state TEA and analyze the flows from the (r, s) = (1,1) to the (r, s) = (1, 2) sectors. 
The thesis is organized as follows: 

In chapter 2, I present an introduction to the basic conformal field theory and define 
the Lee-Yang model, and introduce the necessary methods to be used later. 

In chapter 3 I introduce asymptotic states in defect theories and the notion of the 
transmission matrix. Then I determine the coordinate dependence of defect form 
factors. Ey specifying the boundary form factor axioms we postulate the axioms for 
diagonal defect theories. Using an analogy between defects and standing particles we 
subject our axioms to a consistency check. Then we determine the form factors of 
any bulk operator in terms of the transmission factor and the already calculated bulk 
form factors and outline the procedure to calculate the general solution for operators 
localized on the defects. 

Afterward we apply this technology to determine the defect form factors of the Lee- 
Yang model. Ey calculating the dimension of the operators we can map them to the 
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UV Hilbert space of the model. Finally I introduce a method to derive the different 
boundary conditions via the defects. 

In chapter 4, we start with the definition of the lattice theory. We define the face 
weights and the periodic, defect and boundary raw transfer matrices. We derive the 
functional relations they satisfy. Then we analyze the analytical structure of the 
transfer matrices to turn the functional relations into integral equations, obtaining 
the TEA equations in the three models in the massive and the massless cases. 

We start with the trigonometric/conformal case: First we make correspondence be¬ 
tween the UV Hilbert space in terms of Virasoro modes and the zeros of the transfer 
matrix and the paths. Then we analyze the lattice flow in the parameter ^ and 
describe our findings in the three languages: paths, zeros, modes. We repeat this 
analysis for the periodic and defect cases. We derive the integral TEA equations for 
the massless cases. 

Switching into the massive case, and using the elliptic theta functions for the face 
weights instead of the trigonometric ones, we get the same analytic structure and the 
same paths, zeroes , modes. We derive the massive TEA equations for the ground 
state and obtain the source terms. 

Then in chapter 5, I present the basic results which we obtained |128| for the Luscher 
correction terms to the asymptotic Eethe Ansatz energy of an su{2) particle state in 
finite volume by vacuum polarization effects due to wrapping interactions. 

I present the conclusions in chapter 6. 

Finally the appendix presents how to calculate the finite size correction which origi¬ 
nates from virtual particles propagating around the circle. 

Those Luscher correction terms which were presented in the outline are calculated by 
the weak coupling expansion and shown in the appendix. 
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The material presented in this thesis is based on the papers 1231 El M: 

• Zoltan Bajnok, Omar el Deeb: Form factors in the presence of integrable defects, 
Nuclear Physics B 832, p. 500-519, 2010, arXiv:0909.3200 [hep-th] 

• Zoltan Bajnok, Omar el Deeb, Paul Pearce: Lattice, In Progress 

• Zoltan Bajnok, Omar el Deeb: 6 -loop anomalous dimension of a single impu¬ 
rity operator from AdS/CFT and multiple zeta values, JHEP 1101:054,2011, 
arXiv:1010.5606 [hep-th] 




Chapter 2 


The Lee-Yang Model 


We start by considering a massive scattering theory with n types of particles with 
masses rria- One-particle states are denoted by l^)^, where 6 is the particle rapidity. 
The Lorentz invariant normalization implies that a{(^a\(^b)b = ~ db) 

Asymptotic states of the theory are dehned as tensor prodncts of one-particle states 
and are denoted by \ 0 i, 62 , aia 2 ...an) where 6 i > 62 > ■■■■ > 6 ^ for in states and 
61 < 62 <■■■■< 6 n ioi out states. 

2.1 Conformal Field Theories 

Conformal Field Theories are two-dimensional Euclidean held theories, which possess 
invariance under conformal transformations, including scale-invariance. A general in¬ 
troduction to conformal held theories was presented in [761 [77]. CFTs are a special set 
among the QFT’s in the sense that they represent hxed points under the renormal¬ 
ization group how. In statistical physics, they describe huctuations of critical systems 
in the continuum limit [zg. Conformal invariance highly constrains the behavior of 
the correlation functions, and even the operator content of the theory m- 

Due to scale invariance the energy-momentum tensor T^y{x,y) is constrained to be 
traceless. In complex coordinates z = x + iy and z = x — iy this condition is expressed 
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as T^z = Tzz = 0. Conservation of the energy-momentum tensor means that dzTzz = 
dzTzz = 0. Therefore, being holomorphic and anti-holomorphic quantities T{z) = 
Tzz{z) and T{z) = Tzz{z) are chiral and anti-chiral tensor components respectively. 
T{z) may be expanded into its Laurent-series by summing over its modes around 
2 ; = 0 as 

T{z) = Y,^~"~^Ln 

n 

The algebraic operators satisfy the Virasoro algebra: 

\Lm .hm] (n Tri)Lzi-\-m T 

which is a an extension of the symmetry algebra of the classical conformal group. The 
new element c that appears in this expression is the central charge of the theory. 

The value of c restricts the possible representations of the Virasoro-algebra, the oper¬ 
ator content and the spectrum of the theory. The simplest theories are the minimal 
models |78| . which contain a hnite number of primary helds and possess no additional 
symmetries. For those models the central charge is determined from two coprime 
integers p and q: 

^ _ 6 (p - qf 

^ pq 

Theories with q = p + 1 are unitary theories having no negative-norm states. However 
some non-unitary models are also interesting models to study like the Lee-Yang due 
to their simplicity, and the possibility of generalizing the outcomes of their study to 
physical models of higher complexity. 

The LCs generate the conformal transformations associated with z = x + iy with c 
denoting the central charge of the conformal theory. The same algebra holds with 
associated with z = x — iy, the complex conjugate of z. The L and L commute. Each 
operator family is formed of a primary operator $ and its descendants formed by the 
repeated action of and with negative n. Positive values of n annihilate <h. 

The descendant operators are of the form; 
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0 < ii < .... < ii 


0 < ii < •••• < jj 

The descendant operators will have the levels (/, 1) where / = Yl!n=i ^ ~ X]n=i in- 

The conformal dimensions of the descendant operators will be (A, A) = (A$ + /, A$ + 

1 ). 

A scaling operator $ has the conformal dimension (A$,A<i>) which determines the 
scaling dimension and the euclidean spin by: 


A<j) = A<j) + A<j) 


5 (|) — 


One has to pay attention that not all descendants are independent due to the presence 
of degenerate representations that possess vanishing linear combinations of descendant 
operators. 


2.2 CFTs on the cylinder 


The conformal transformations 


= exp ~ 


2 ; = exp 
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can be used to map the complex plane into a cylinder of spatial circumference L , 
where x and r are the spatial and the imaginary time coordinates, respectively. With 
this transformation the Hamiltonian operator will be defined as: 

27r / — c \ 

H= dx Trrix^T = 0) =—{^Lo + Lo -—j (2.4) 

In minimal models the Hilbert-space is given by 

^ = (2.5) 

h 

where V and V denote the irreducible representation of the left and right Virasoro 
algebras with highest weight h. 

2.3 Perturbing CFT’s 

Conformal field theories represent statistical physical or quantum systems at critical¬ 
ity. However, they can be also used to approach noncritical models. 

If there is only one perturbation present, which only brakes a subset of the conformal 
symmetries, the theory may still possess an infinite number of conservation laws, and 
it may remain integrable HSllHnilEI]. We assume that the theory defined by the action 

A = Acft d J <h(a;) (2.6) 

is integrable. Scale-invariance is broken by the perturbation. This action can define 
a massless or a massive perturbation depending on the original CFT and the nature 
of the perturbation. 

The correspondence between perturbed CFTs and scattering theories can be in¬ 
spected by several ways. 

• Zamolodchikov’s counting argument ESI. 


CHAPTER 2. THE LEE-YANG MODEL 


16 


• Thermodynamic Bethe Ansatz [H21IS3], which determines the finite size depen¬ 
dence of the vacuum energy in terms of the S-matrix of the theory. It allows to 
predict the central charge of the CFT and the scaling dimension of the perturb¬ 
ing field. 

• Truncated Conformal Space Approach (TCSA) which can be used to numeri¬ 
cally determine the low-lying energy levels of the finite size spectrum. One can 
identify multi-particle states and test the phase shifts Sab{0)- 


2.4 The Lee-Yang model 

The non-unitary minimal model M 2,5 has central charge c = —22/b and a unique 
nontrivial primary field $ with scaling weights A = A = —1/5. The field $ is 
normalized so that it has the following operator product expansion: 

d)(z, z)d)(0,0) = 0) + {zzfH + ... (2.7) 

where I is the identity operator and the only nontrivial structure constant is 

C = 1.911312699- X i 

The Hilbert space of the conformal model is given by 

Uly = ® Vh 

h=o-1/5 

where Vh (V/x) denotes the irreducible representation of the left (right) Virasoro algebra 
with highest weight h. 

The off-critical Lee-Yang model is defined by the Hamiltonian 

H^^^ = Hcft + iX [ dx^{ 0 ,x) (2.8) 

Jo 


CHAPTER 2. THE LEE-YANG MODEL 


17 


where He ft is the conformal Hamiltonian. 

When A > 0 the theory above has a single particle in its spectrum with mass m. 
The S'-matrix reads [SI] 


^ ^ sinh 9 + i sin ^ 
SLYiO) = 3 


sinh 9 — i sin ^ 


(2.9) 


and the particle occurs as a bound state of itself at 6^ = 2 m/II with the three-particle 
coupling given by 


= - 2^3 


where the negative sign is due to the non-unitarity of the model. 



Chapter 3 


Form Factors in Presence of Defects 

3.1 Defect form factors 

In this section we present the axioms for the matrix elements of local operators be¬ 
tween asymptotic states. To shorten the discussion we introduce Zamolodchikov- 
Faddeev (ZF) operators in order to describe both the multiparticle transmission pro¬ 
cess as well as the properties of the form factors. 

3.1.1 Asymptotic states and transmission matrix 

Multi-particle asymptotic states in integrable bulk theories can be formulated in terms 
of the ZF operators as 


\ei,...,en) = A+{e,)...A+{en)\o) 

All particles have different momenta pi = msinh. 6i, thus in the remote past they are 
not interacting and form an initial state 6 i < ■ ■ ■ < 6^. When time evolves they 
approach each other and after the consecutive scatterings they rearrange themselves 
into the opposite order; 
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i<j 

Here S{9) is the two particle scattering matrix which satisfies unitarity and crossing 
symmetry 


S{-9) = S-\9) ; S{i7i -9) = S{9) 

This mnlti-particle scattering process is easily formnlated with the ZF algebra: 


A+{9,)A+{92) = S{9i - 92)A+{92)A+{9,) + 27r6{9, -92- tn) (3.1) 

where we extended their definition for imaginary 9 by postnlating the crossing prop¬ 
erty 

A{9) = A+{9 + in) (3.2) 

Once defects are introduced we have to make a distinction whether the particle arrives 
from the left (H) or from the right (B) to the defect. These particles can be even 
different from each other as they live in different subsystems. A multiparticle state is 
then described by 

|0i,..., ...,9^)= A+{9 ,)... A+{9n)D+B+{9n+i) • • • H+(0™) |0) 


where the ZF operators B^ create particles on the right of the defect and satisfy 


similar defining relations to (3.1) with possibly different scattering matrix. Yet, for 
simplicity, we restrict our discussion to the case when the two subsystems are identical 
with the same scattering matrix. Observe however, that this does not imply space 
parity invariance, since the defect may break it. In the initial state rapidities are 
ordered as > ■ ■ ■ > > 0 > 9n+i > ■ ■ ■ > 9m- The final state, in which all 

scatterings and transmissions are already terminated, can be expressed in terms of 
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the initial state via the multiparticle transmission matrix. 


n 


m 


\e,,...,e^-e^+^,...,em) = \{s{e,-e^)\[T_{e,) J] t +(- 0,)|0 



i<j 




Due to integrability it factorizes into pairwise scatterings and individual transmissions: 
T-{9) and T+(—0). We parametrize T+ in such a way that for its physical domain 
{9 < 0) its argument is always positive. Transmission factors satisfy unitarity and 
defect crossing symmetry |34( 


r+(-9) = r-‘(9) i r_(9) = r+(»,r - 9) 


(3.3) 


The multiparticle transition amplitude can be derived by introducing the defect op¬ 
erator and the following relations in the ZF algebra: 


yl+(0)D+ = T_(0)D+S+(0) ; D^B+{-9) = T+(0)Zl+(-0)D+ 


A defect is parity symmetric if T_{9) = T^{9). Clearly A+(0 = 0) satisfies the 
properties of with T-{9) = S{9) = T^{9) . Thus a standing particle can be 
considered as the prototype of a parity symmetric defect. 

3.1.2 Coordinate dependence of the form factors 

The form factor of a local operator 0{x,t) is its matrix element between asymptotic 
states: 



where the adjoint state is defined to be 


^ ^ • • • B{9l,^^)DA{9\) ... A(0;) 
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Strictly speaking the form factor is defined only for initial/final states (i.e. for de- 
creasingly/increasingly ordered arguments) but using the ZF algebra we can generalize 
them for any values and orders of the rapidities. 

The multiparticle asymptotic states are eigenstates of the conserved energy. This fact 
can be formulated in the language of the ZF algebra as 

[H, kl+ (6)] =m cosh 6 Zl+ (0) ; [H, D+] = 

In the second equation we supposed that the vacuum containing the defect has energy 
cd- Classical considerations together with the topological nature of the defect suggest 
the existence of a conserved momentum with properties 

[P,A+{e)]=msmheA+{e) ; [P,D+]=pdD+ 

Thus, opposed to a general boundary theory, the defect breaks translation invariance 
by having a nonzero momentum eigenvalue pr, and not by destroying the existence of 
the momentum itself. As a consequence the time and space dependence of the form 
factor can be obtained as 


• • • ,C'+i;Co • • • ... ,9n;0n+i,... ,9m) = 

^itAE—ixAP (nf n' . n' 0^ \0 0-0 0 \ 

^ in' ^rn''){n^rn) * * * ’ ^n/+l ’ ^ nl'i • • • ? ”n+l ? • • • ? ^n+m) 


where NE = m(X)jCosh6'j — Xly cosh6''/) and AP = m(X)jSinh6'j — sinh6*'/). 

The very same simple space and time dependence can be seen in a theory without 
the defect and it is substantially different from what we would expect from a general 
boundary theory where even the one point function has a nontrivial space-dependence. 
These considerations remain valid for operators inserted at the defect 0{t) = O{0,t), 


too. 
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3.1.3 Crossing transformation of defect form factors 

The properties and analytical structure of the form factor F(n'm'),{nm) can be derived 
via the reduction formula from the correlation functions similarly to the boundary 
case [29]. Instead of going to the details of the calculation of I2SI we note that all 
equations follow from the dehning relations of the ZF algebra and the locality of the 
operator [0(0, 0), d+(6*)] = 0 except the crossing relation. It reads as 


T^o ml 

^ {nl,ml){n,m)Tnl+mn • 




nl+mn 


) ^nl+11 ^nn • ••) 1 ^ 1 ) • • • ) ^n+li ^n+m) 

■■■1 ^n/+l) ^nli •••) 1^1) • • • ) ^nj ^n+li ^n+m—l) 


and can be obtained as follows: We fold the system to a boundary one: B^{9) -H- 
B^{—9), and consider A'^ and B^ as creation operators of two different type of 
particles which scatter trivially on each other. Now we apply the crossing equation of 
B^ for the resulting boundary form factor. If we fold back the system to the original 
defect theory we obtain the defect crossing equation above. 

By analyzing the crossing equation of the particle A'^ instead of B^ we obtain 


^{nl,ml){n,m) (.^nl+mli ^nl+1 > ^nli • ■ ■ 5 1^1; ■ • • 5 ^n+1 1 ^n+m) 

^{ni+l,ml){n-l,m)i^nl+mn ^nl+li ^nli ~ • 5 ^n+1) ■■•5 ^n+m) 


This crossing equation can also be obtained from (3.2). Using any of the crossing 


equations above we can express all form factors in terms of the one-sided form factors: 

• • • 1 ^ni ^n+1) ■■•5 ^n+m) •= ^{0,0){n,m){i • ■ ■ ) ^n+lj ^n+m) 


The properties of this form factor follows from the ZF algebra relations and from the 
crossing relations and we postulate them in the next subsection as axioms. 
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3.1.4 Defect form factor axioms 

The matrix elements of local operators satisfy the following axioms; 

I. Transmission; 



By means of this axiom we can express every form factor in terms of the elementary 
one 

It satishes the further axioms; 

II. Permutation; 



III. Periodicity; 


,.. .9n,... ,0i - 2in) 


Ff(9i,92,...,9„) = jy(92 
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The physical singularities can be formulated as follows. 

IV. Kinematic singularity: 



V. Dynamical bulk singularity: 


-iRes0i=0F^_^_2{9' + —,9 — 


.,9r^) = TF^^,{9,9,,...,9r,) 




where T is the 3 particle on-shell coupling. 

VI. Dynamical defect singularity: 


tRes0=,uF^^,{9,, ...,9n,9) = igF^{9,, ...,9^) 
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where g is the defect bound-state coupling. 

A few remarks are in order; Although the form factor axioms (II-V) are the same as 
the axioms of the form factors in a theory without the defect, the axioms (I,VI) are 
different and in general we will have different solutions. An exception is the invisible 
defect T_(0) = 1 when we recover the usual form factor equation providing a consis¬ 
tency check for our axioms. Another consistency check can be obtained by considering 
a standing particle as the defect. Then T±{6) = S{6) and the two additional axioms 
become part of the old ones: (I,VI) will be special cases of (II,V), respectively. 

3.2 Form factor solutions, two point functions 

In this section we determine the solutions of the form factor equations for operators 
localized in the bulk and at the defect. For operators localized in the bulk the solutions 
can be built up form the bulk form factors and from the transmission factors. Using 
the defect form factor solution we determine the spectral representation of the two 
point function for the situations when the operators are localized on the same or on 
the opposite sides of the defect. Finally, for operators localized on the defect we 
outline the strategy for the general solution. 

3.2.1 Bulk operators 

The form factor axioms for are the same as in the bulk so we expect to use the bulk 
form factor solutions. Clearly we have to make a distinction whether the operator are 
localized on the left, or on the right of the defect. If the operator is localized on the 
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left then particles arriving from the left can reach the operator without crossing the 
defect. Since the defect is topological changing its location does not alter the form 
factor (as far as we do not cross the insertion point of the operator). Shifting the 
defect faraway we expect to obtain the form factors of the bulk theory. 



Thus we can conclude that for the initial state 6i > ■ ■ ■ > 6n > 0 the defect form 
factor coincides with the bulk form factor. Let us denote the solutions of the bulk 
form factor equations by Bn{0i ,..., 6n)- Then we claim that for an operator localized 
on the left for short) we have 


F<{0u...,»n) = B„{eu ■■■,»„) = 




(3,4) 


By using the transmission axiom and the crossing relation we can express all other 
matrix elements in terms of the bulk matrix element and the transmission factor. If the 
operator is localized on the right of the defect {O^) then, by similar argumentation, 
we expect the defect form factor to coincide with the bulk form factor for particles 
coming from the right. Those initial states have the ordering 0 > > • ■ • > and 

the form factor is then 


The solution for the elementary defect form factor for operators localized on the right 
thus turns out to be 


... , 9 „) = nr-(9i)B„(9i, ..., 9„) 


(3,5) 
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which satishes all the bulk form factor axioms but does not coincide with the bulk 
form factor solution. Having calculated the form factor solutions we use them to 
construct the two point functions, which for operators localized on the opposite side 
of the defect will be intrinsically different from the one without the defect. 

We analyze the following two point function 

{]\Oi{Xi^ ^ 2 )!; ) 

where we denote by |;) the vacuum of the defect theory. Formally |;) = F)''‘|0). Now 
we insert the resolution of the identity. It can be composed both from initial and from 
hnal states and for dehniteness we choose initial states. It is instructive to list the 
possible states. If we have no particles we have only the vacuum; |;). One particle 
states can be of two types, depending on whether the particle arrives from the left or 
from the right: |0;) for 0 > 0 and |; 9) for 0 < 0. A general N = n + m particle state 
| 6 'i,..., On] On+i ,..., On+m) with 6*1 > ■ ■ ■ > > 0 > On +1 > '' • > 6',^ has to cover all 
possible cases ranging from n = 0 to n = iV . The two point function then can be 
written formally as 

OO 

v=o 

We have to specify the integration ranges for the multiparticle state N. Originally 
we have to integrate only for the multiparticle momentum range of the initial states. 
If we exchange the order for di > 62 to the nonphysical 62 < 0\ then the form factor 
of 0\ picks up a factor SiO^ — 62 ) while that of O* the inverse factor S*{ 6 i — O 2 ), so 
the integrand is a symmetric function. For each integration with > 0 we also have 
an analogous integration for —0i < 0 . Their contributions differ by a factor T+(6'i) 
for the form factor of Oi and by the inverse T^(Oi) for O 2 . As a consequence we can 
express the correlator in terms of the elementary form factors Fn as: 
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(; \0i{xi,ti)02{x2,t2)\;) = 

■spoo TOO roo d^j^Oi(a a a '|*pi_E(n)(ti-t2)-iP(n)(xi-a;2) 

2-^n=0 n\ J—CO 27 r * * * J—co 2 'k h V I 5 * * * ? nj^n \”1? • ■ • ? n) 

Although we transported the operators Oi and O 2 into the origin, the form factor 
solutions remember whether the operators are localized on the left or on the right of 
the defect. 

If both operators are localized on the left, {xi < 0, 0:2 < 0) then the elementary form 


factors are the same as the bulk form factors (3.4) and we can conclude that the two 


point function is exactly the same as the bulk two point function 

(; |C>l(xi,tl)C>2(x2,t2)|;) = 

J_ roo ^ roo d^ToOi(f) Q ft \* iE{n)(ti-t2)-iP{n){xi-X2) 

l^n=0n\J-oo 2 -k ' ' ' J-oo 2 -k , • • • , , . . . , t7„j e 

This is intuitively clear: we can transport the defect to infinity without crossing any of 
the insertion points thus leaving invariant the two point function. When the defect is 
at infinity it does not influences the two point function which then has to be the same 
as in the bulk. The same result can be obtained when both operators are localized 
on the right of the defect. 

If the operators are localized on different sides of the defect (xi < 0, 0:2 > 0) then 


additionally to (3.4) we also have to use (3.5). As a result the two point function is 


expressed in terms of the bulk form factors and the transmission matrix T_ (6*) as 


Er=o h fZo ■■■fZo ^TA0n)BOReu ..., ..., 

This is the main result of this section. This formula shows how the correlation function 
can be calculated in the presence of an integrable defect in terms of the transmission 
factor and the bulk form factors. It can be generalized to any correlators localized 
in the bulk using the resolution of the identity together with the exact form factors 


(3.4) and (3.5). It cannot be applied, however, for operators localized on the defect. 


which is the subject of the next subsection. 
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3.2.2 Defect operators 

We have seen that although the minimal form factors Fn are subject to the same 
requirement as the bulk form factors they are not necessarily the same. In this 
subsection we develop a general methodology to determine the defect form factors. 
Let us analyze them for increasing particle numbers: 

The hrst form factor is the vacuum expectation value of a defect held 

{■,\0{t)\-) = Fo 

The one particle form factor is dehned to be 

{■\o{t)\e-) = F,{e) ; {■\o{t)\-e) = T_{e)F,{e) 

Contrary to the bulk case it has a nontrivial rapidity dependence: it is not natural to 
take Fi to be a constant. In a parity invariant theory for a parity symmetric operators, 
for example, we have Fi{6) = T_{—6)Fi{—d). If parity is broken then Fi{6) can be 
an arbitrary defect condition-dependent 27rf-periodic function. The only restriction 
came from the defect bound-state axiom (VI): it must have a pole at iv whenever 
T_{6) has a pole corresponding to a bound-state. Let us denote the minimal function 
which satishes this requirement by d{6). The general form of the one particle form 
factor is then 

{■,\0(m} = d(ff)Q°(x) ; x^e" 

where d{6) depends on the defect condition, while Q{x) depends on the operator we 
are dealing with. 

The two particle form factor must also have a singularity at iu and additionally it 
satishes the bulk form factor axioms so we expect it to be written into the form 


.^2(^1, ^2) — d { 9 i ) d { 92 ) —6^2) 
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where fmin{(^) is the minimal solution of the bulk two particle form factor equations 

= S{e)fmin{-0) ] - ^) = + 0) 

Taking into account the general parametrization of the bulk and boundary form factors 
together with the dynamical and kinematical singularity axioms we parametrize our 
minimal defect form factors as 

Fn{Oi,...,0n) = W d{9i) _ _ _ ^ 3 ;^) . Xi = e^^ 

. ^ . Xi ~r Xj 

^ i<j 

where Q{xi,... ,Xn) is a symmetric function expected to be a polynomial, if there 
is no bulk dynamical singularity. If there is such a singularity we have to include 
the corresponding singularity into fmin ■ The dependence on the defect condition is 
contained in d{6) with possible defect bound-state singularities, while the dependence 
on the operator is in Q. If for instance the defect is the invisible defect with T± = 1 
then d = 1 and we recover the solution of the bulk form factor equation as it should 
be. From the kinematical singularity equations recursion relations can be obtained 
among the polynomials Qn +2 and Qn- 

3.3 Model studies 

In this section we analyze the solutions of the defect form factor axioms for the free 
boson and for the Lee-Yang models. 

3.3.1 Free boson 

The purely transmitting free bosonic theory was analyzed in EH as the limiting case 
of the sinh-Gordon theory. The Lagrangian of the model reads as 
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C = Q{-x) 
6{x) 


m 




+ 0(x) 


1/0 ^ \2 
-(d^^+y - 


2' ^ ' 2 

— $_<i>+ + m [(coshju) + 2(sinhju) $+$_] 


where are the helds living on the right/left part of the defect and p is a free 
parameter. By varying the action we obtain the free equation of motion in the bulk 


(d^ - (92)$± = 


and the defect conditions; 


±(9t<h± =F = m(sinh ju) <l>± + m(cosh ju) <l>=p 


Since are free helds, they have an expansion in terms of plane waves and cre¬ 
ation/annihilation operators 


^±(x,t) = 


p dk 1 

/_oo 27r 2u{k) 


+ ; uj{k) = Vk^ + 


where the a, operators are adjoint of each other with commutators: 


[a±{k), a^{k )] = 27r2u{k)S{k — k ) 


They are not independent, the defect condition connects them as 


a±(±/c) = Tzp(/c)a=p(±/c) 


T^{k) 


msinh/i =f iuj{k) 
m cosh fi — ik 


k>0 


This shows that the defect is purely transmitting, that is we do not have any rehected 
wave. The transmission factor in the rapidity parametrization [k = mcisinh6') can be 
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written also in the following form; 

Tje) = - f + 2 ) _ I + wx 

sinh(| + f + f) l-wx 

where w = and x = e^. Sometimes we also use w = w~^ and x = x~^. In the next 
subsection we will set m = 1 and use dimensionless quantities. 

3.3.1.1 Form factors 

In the Free boson model, we have the advantage that we can explicitly calculate the 
form factors of all the operators, and then check that they satisfy the defect form 
factor axioms. Additionally we can also conhrm that we have as many polynomial 
solutions of the axioms as many local operators exist in the theory. 

We work with the Euclidean version of the theory [t = iy) and introduce complex 
coordinates z = y + ix,z = y — ix. We use the explicit expressions of $±( 2 :, z) above 
to calculate the form factors. The one particle form factors turn out to be: 

Ff- = (O|$_(;2,^)|a_(0)) = 

Ff+ = (O|$+(;2,f)|a_(0)) = e“^T_(0) 

from which it is easy to calculate the defect form factors of the derivative of the 
elementary helds; 

(O|9"d>_(O)|a_(0)) = x^ ; (O|a"$_(O)|a_(0)) = x^ 

{O|a’"$+(O)|a_(0)) = x^T_{e) ; (O|(9’^d>+(O)|a_(0)) = x^T_{e) 

We can unify the notation by (9“" = d'^. It is instructive to see how we can re¬ 
cover these form factors from the solution of the form factor axioms. Now using the 
parametrization of the form factors in terms of d{d) , we know that at level 1 the 
solutions of the form factor axioms have the form; 



CHAPTER 3. FORM FACTORS IN PRESENCE OF DEFECTS 


33 


Thus if we choose 


we obtain 


Fi(0) = d{9).Q,{9) 


d{9) 


1 

1 — wx 


Qi^ (9) = x'^{l — wx) ; (9) = x^{l + wx) 

Naively it seems we have less polynomial solutions of the form factor equations as 
operators; We have extra relations among the form factors originating from 

= ; (9$_ - (9$+=«;($+ +$_) ; (9<h_ + 9$+=«;($+-$_) 

However, these relations are satished due to the bulk equation of motion and the 
defect conditions. Note that the form factor solutions are even more simple in terms 
of 0 = <h+ + <h_ and 0 = Actually 0 is not independent since 0 = wd(j). 

Their form factors read as: 

QT^ = X- ; QT^ = wx^^^ 

The general multiparticle form factor as calculated from the explicit solution of the 
model reads as 

(0| : _: \9i,..., 9^) 

= x^" ... x^T{xk+i) ■. .T{xn) + permutations 

where Xi = e^\ In the form factor bootstrap the general parametrization is 

N 

Fn{'i 6 ^ 1 , • • •, 6 * at ) = d{9i)Q]\[{xi ,..., xm) 

i=\ 
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Thus we can read off the corresponding form factor solution 

Qn = (1 — wxi )... (1 — wxk){l + wxk+i )... (1 + wxn) + permutations 

Since the scattering matrix in the free boson theory is trivial S' = 1, the form factors 
of different levels are not connected to each other, and in this way we solved the 
theory completely. 

In terms of the field (j) the form factor solutions are exactly the same as in the bulk 
free bosonic theory: 

Qn = ... x^ + permutations 

thus we have exactly the same number of solution of the form factor axioms as many 
independent local operators in the theory. 


3.3.2 Defect scaling Lee-Yang model 

The scaling Lee-Yang model can be defined as a perturbation of the At ( 2 , 5 ) conformal 
minimal model with central charge c = — ^. It contains two chiral representations 
of the Virasoro algebra, Vq, Y with highest weights 0 and — respectively. The 
fusion rules can be summarized as: = A*q = 1 and Nl^ = 1 for z = 0,1 and all 

others are zero. The Hilbert space on the torus corresponds to the (diagonal) modular 
invariant partition function and contains modules corresponding to the Id and the 
<h( 2 :,z) primary fields with weights ( 0 , 0 ) and (—|, —|): 

Ti = Vq ®Vq -\- Vi (3-6) 


The only relevant perturbation by the field $ results in the simplest scattering theory 
with one neutral particle of mass m and scattering matrix | 12 ] 


s{e) 


sinh 9 T i sin | 
sinh 9 — i sin | 
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The pole at 0 = y (with residue T^) shows that the particle can form a bound-state. 
The relation 

however, implies that the bound-state is the original particle itself and the bulk boot¬ 
strap is closed. 


3.3.3 Integrable defects 


Two types of topological defects can be introduced in the At ( 2 , 5 ) minimal model. They 


can be considered as operators acting on the bulk Hilbert space (3.6) commuting with 


the action of the left and right Virasoro generators. They have to act diagonally on 


each factor in (3.6) and satisfy a Cardy type condition. This leads to two choices which 


can be labeled by the same way as the bulk fields: (0,0) and (1,1). After making a 
modular transformation the defect is inserted in space and the corresponding Hilbert 
space can be described as 






For the topological defect labeled by (0,0) the Hilbert space turns out to be 


^(0®) = Ho 0 Ho + Hi 0 Hi 


and coincides with the bulk Hilbert space. This defect is the trivial (invisible) defect. 
For the other defect labeled by (1,1) we obtain 

rAi) = + VI ® Ho + Ho 0 Hi + 2 Hi 0 Hi 


For each of the representation spaces we associate a primary field Id, ^) 

and ^+(z,z) with highest weights (0, 0), (—1/5,0), (0,—1/5), (—1/5,—1/5), respec- 
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tively. The non-chiral fields ^±{z,z) can be considered as the left/right limits of the 
bulk field $ on the defect. 

The bulk perturbation by $ in the defect conformal field theory does not break inte- 
grability. In the case of the trivial defect the transmission factor is simply the identity 
T = 1. In the case of the defect labeled by (1,1) we can introduce a one parameter 
family of defect perturbations by properly harmonizing the coefficients of the (p{z), 
(p(z) and ^{z,z) terms. We plan to analyze this issue in a forthcoming publication. 
Related investigations with only defect perturbations can be found in jl^. The re¬ 
sulting theory is integrable and can be solved by exploiting how the defect acts on 
integrable boundaries |3Z]. In the calculation the bootstrap relation 


?7r ?7r 


(3,7) 


was used together with defect unitarity and defect crossing symmetry (3.3) to fix the 
transmission factor as 


r_(()) = |4+i]|6-i] 



siiih(| + ) 

sinh(| + - *f) 


(3,8) 


Actually the inverse of the solution is also a solution but the two are related by the 
b ^ 6 + b transformation. 

We also note that the defect with parameter 6 = 3 behaves as a standing particle 
both from the energy and from the scattering point of view. 


3.3.4 Defect form factors 


In this subsection we apply the general method developed in Section 3 to determine 
the form factors of the defect Lee-Yang model. The form factor can be written as 


Fn{eu...,en) = Hnl[d{e,)l[ 


Anin(^* 






i<j 


Xi -h Xj 
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The minimal solution of the two particle form factor equation is well-known reads as 

lig: 

/miri(®) = >'(“ - 9) + 0 ) 

where 

f dx sinh f sinh f sinh I 

v{6) = exp < 2 / —e - -^^ 

{ Jo X sinh X 

We also included the pole corresponding to the dynamical singularity equation by the 

denominator. We choose the normalization of the form factors as in the bulk 

^ f 3a y 
" 4v/3 l^ 2 iu( 0 ) J 

Qn{xi,, Xn) is expected to be a symmetric polynomial in Xi and Xi. 

Let us turn to the determination of d{6). Due to the defect dynamical singularity for 
Fn{0i,..., On) the defect dependent term d{9) must have a pole whenever T-{0) has a 
pole. Similar equation is valid for 9i,..., 9n) 

defect bound-states poles of Tj^{9) . We will take into account that the transformation 
606 — 6 exchanges T_{6) with T^{9) and we expect that it acts in a similar way on 
the form factors. The minimal solution with these requirements turns out to be: 


di9) = 


4sinh(2-f 3 ^(6 — 5)) sinh(|-f (6 — 7)) + 2 cos{^ — i9) \/3 + xu + x 


where we introduced z/ = and V = z/“^. This function satisfies two relevant 
relations: 

d{9 + in)d{9) = 


and 


1 — 2 cos(y — 2 i 6 *) 1 — 


ITT ITT 1 1 

d{9 + -)d{9--) = 


3 2cos(^—-Z0) XV P X V 


Singularity axioms generate recursive relations between the polynomials. The kine- 
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matical recursion relation is given by: 


Qn+2( ^ 1 ^ 1 • •• 1 ^n) Dn{x ^ ^ Xn)Qn{A^'i '"i ^n) 


with 


Z)„(x,xi, ...,Xn) = + n 

/ -V. 

1 


2x{oj — uj) 


J^(xa; + Xja;)(xa; — XiUj) — J^(xa; — Xiw)(xa; + Xjw) 


^i=l 


2=1 


where we introduced u = , u = u while the the bound state recursion relation 

is ; 


n—1 

Qn+l{xU,XU},Xi,...,Xn-l) = (z/X^ + Z/) (x + Xi)(5„(x, Xi, ..., X„_l) 

2=1 


Now I proceed to solve these recursions. 


3.3.4.1 Solutions 

Since Q„(xi,x„) is a symmetric polynomial, it is useful to introduce the elementary 
symmetric polynomials a^^\xi, ...,Xn) dehned through the generating function; 

n n 

JJ(x + Xi) = ...,Xn) 

i=l k=0 

By means of these functions the kinemetical recursive relation for Qn reads as: 

(-1)”+^Q„+2(-X, X, Xi, ...,Xn) = - X^ + V~‘^)Dn{x, Xi, ...Xn)Qn{xi., X^) 


n k 

£>„(x,Xi, ...Xn) = ^ ^ 

k=l m=l,odd 


2(n-k)+m (n) (n) / -\\k+l 

sinCf) ^ ’ 


We are going to hnd the form factors of the operators *h±(x, z), <p(x), ip{z) and their 
descendants. We can choose <h±as the defect limits of the right/left bulk helds, thus 
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we know already all of their form factors. Taking into account the explicit form of 
d{9) together with T_(6*) we hnd 

Qf~ = I'd + Pa + Vs ; = va + Pa — Vs 

For two particle form factors we get 

Q2~ = ai{v‘^a2 + VSvai + aiai + 1 + VSP&i + p‘^a2) 

Q2^ = ai(v^a2 — VSvai + ciiai + 1 — VSPai + p'^a2) 

where we used the solution of the bulk form-factor equation Qf = ai. They satisfy 
the dynamical recursion relations. The asymptotics of the solutions for x —)■ ±cx) 
reflect the dimensions of the helds (—|,—|). We would like to describe two more 
chiral helds (f(z) and (f(z) with dimensions (—1,0) and (0,—|). The corresponding 
solutions at level one turn out to be 


Qf = ai ; Qf = ai 

They are related by the x -H- x~^ transformation. Using our recursion relations we 
hnd the related solutions at level 2 

Q^ = ai(vaf + P) ; Q^ = ai(Paf + v) 

We summarize the FF solutions of the primary helds up to level 2 in this table: 


Operator 

Qi 

Q 2 


lyai + Pai + VS 

ai(v^a2 + VSvai + aiai -|- 1 -|- VSPai + P^a2) 

$+ 

uai -|- Pai — Vs 

ai(v^a2 — VSvai + aiai -|- 1 — VSPai + p‘^a2) 

F 

0-1 

ai{vaf + P) 

F 


ai{Pal + v) 
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Together with the Identity operator, we have 5 independent operators in this model, 
and with their descendants they describe the full spectrum of the solutions. How¬ 
ever, the one-to-one matching remains an open problem due to the large number of 
operators and their descendants in this model. 

We also list the first and second order descendant operators and their corresponding 
solutions at levels 1 and 2 below; 


• First Order Descendants: 


Operator 

Qi 

Q 2 

d(j)+ 

ai{vai + uai — ^/3) 

af(v^a2 — y/Svai + cricri -|- 1 — \/3v~^d'i + v~^d'2) 

d(j)+ 

ai{vai -|- Pai — ^/3) 

di(Ti(n^(j2 — \/3vai -|- aidi -|- 1 — \/3v~^ai + v~^a2) 

(90_ 

(Ti {vai + z/cfi -h a/3) 

aj{v‘^a2 + a/Sucti -h aidi + 1 + a/Sw^^i -F ^“^^ 2 ) 

(90_ 

ai{vai + uai + a/3) 

di(Ji(n^(T2 -|- VSvai -|- ai&i -|- 1 -|- \/3v~^ai P v~‘^a2) 

dip 


crf(vaf + v~^) 

Bp 


aid-i(v~^d-f + v) 


• Second Order Descendants: 


Operator 

Qi 

Q 2 

d‘^B+ 

afOai + v~^ai — a/S) 

(jf(w^CT2 — a/SuCTi -I- CTifTi -I- 1 — \/3v~^ai + V“^CT2) 

dd(j>+ 

aiaiRai + v~^ai — a/3) 

(71(Tj(w^CT2 — a/3W(Ti -I- (Ti(7i -|- 1 — A/3w“^(tl -|- ^“^^ 2 ) 

d'^(j)+ 

aKvai + v~^ai — a/3) 

o-^fTi(w^CT2 — VSvai + -I- 1 — A/3v“^(ti -|- v~^a2) 


crj(wi +v~^ai + a/3) 

af{v^a2 -t a/Swcti -I- ai&i -I- 1 -I- \/3v~^ai + w“^CT2) 

dd(j>- 

aiuRvui + v~^ai + a/3) 

aia‘f{v‘^a2 -t a/3vcti -I- cjiai -I- 1 -I- V5v~^ai + v~'^aA 


Cilwi +v~^ai + a/3) 

a\a\{pP'U2 -t a/Svcti -I- a\a\ -I- 1 -t V5v~^ai + v~^a2) 



(Tiivaf -1- a“^) 



aial{v~^al + v) 

T 

0 

ai{ai + a2){v‘^a2 + a/3wi -I- aiai -I- 1 -I- Civ~^ai -I- v~‘^a2) 
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3.3.4.2 Parity Symmetry 

In this part we analyze how the parity transformation acts on the form factor solutions. 
The action of the parity operator P on the operators can be written as 

POp-i = 


The action on the form factors is 

P{0|C'(0)|9; ) = {0|C''’| ;-9) 

Checking for <h_and <h_|_up to level 2, we hnd that they are parity even with 

P{0|4_(0)|«i ) = {0|4+(0)| ;-«) 
while on the contrary, (p and p are parity odd with 

P(0M0)|9; ) =-(0^0)1;-9) 

To conhrm these parity properties of the primary helds one should work out the defect 
Lee-Yang conformal held theory. 

3.4 Boundary form factors via defects 

In this section we intend to illustrate how defects can be used to generate new bound¬ 
ary form factor solutions from old ones. The underlying fusing idea for the rehection 
matrices can be explained as follows: Suppose we place an integrable defect with 
transmission factor T_{6) in front of an integrable boundary with rehection factor 
R{0), which satishes unitarity and boundary crossing unitarity; 

771 ITT 

R{-e) = R-\e) ; R^--e) = s{2e)R{- + e) 
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If we fuse the defect to the boundary the composite boundary system will be integrable 
and will have reflection factor 


R{9) = T+{9)R{9)T_{9) 


which, due to the defect unitarity and crossing equations, will satisfy boundary uni- 


tarity and crossing unitarity. This idea has been used to calculate the transmission 


factors from the already determined reflection factors R, R in the sinh-Gordon and 
Lee-Yang models in EH. In contrast, here we would like to use the fusion idea to 
generate new form factor solutions from old ones. For this purpose we suppose that 
we determined already the boundary form factors F^{9i ,..., 9n) of a boundary op¬ 
erator O. It satishes besides some singularity axioms the following requirements: 
permutation 



reflection 



and crossing reflection 



We claim that the fused form factor 


n 


, 9 „) = n T-mF°{eu ..., 9 „) 


(3.9) 


satishes the boundary form factor axioms of the fused boundary corresponding to 
the rehection factor R. Let us analyze them one by one. Since the extra factor is 
symmetric in 9i the permutation axiom is trivially satished. To show the rehection 
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property we use defect unitarity 

R{e) = T+{e)R{e)T_{e) = T_{-e)-^R{e)T_{e) ( 3 . 10 ) 

while for the crossing reflection we use defect crossing symmetry: 

R{i7r -e) = T+iiir - e)R{i7r - 0)T_(z7r -6)= T_{e)R{i7r - e)T_{2m - 9)-^ (3.11) 


Now multiplying both sides of the reflection and crossing reflection equation by 


])([^T_(6'j) and using (3.10) and (3.11) the claim follows. Similarly one can show 
the satisfaction of the singularity axioms. 

By this method form factor solution of a given boundary can be used to generate 
form factor solutions for the fused boundary. It is practically useful if we can follow 
the identification of the operators under the fusion procedure. This is the case for 
example if the operator in the UV limit commutes with the defect. Say for example if 
in the Lee-Yang model we determine the form factors of the operators of the identity 
module on the trivial boundary, then by the fusion procedure we can generate the 
form factors of the same module on the fused (j) boundary. 




Chapter 4 


Spectrum of Lee-Yang model in finite 
volume 


We consider the non-unitary Yang-Lee minimal model Ai{2,5). It is obtained 1^2] 
as the continuum scaling limit of the A 4 lattice model of Forrester-Baxter BHI in 
Regime III with crossing parameter A = ^. We consider it in three different hnite ge¬ 
ometries: on the strip with integrable boundary conditions labeled by the Kac labels 
(r, s) = ( 1 , 1 );( 1 ; 2 ), on cylindrical geometry with either periodic boundary condi¬ 
tion or by including an integrable purely transmitting defect. We then apply </ 9 i 3 
integrable perturbations both on the strip’s boundary and the defect and describe 
the flow of the spectrum. Introducing an additionally $ 1^3 integrable perturbation in 
the bulk we can go off-critical and determine the hnite size spectrum of the massive 
scattering theory in the three geometries, via thermodynamic Bethe ansatz (TBA) 
equations. We derive these equations for all excitations by solving, in the continuum 
scaling limit, the TBA functional equation satished by the transfer matrices of the 
associated A 4 lattice model of Forrester and Baxter in Regime III. The excitations are 
classihed in terms of simple (m, n) systems. The excited state TBA equations agree 
with the previously conjectured equations in the strip and periodic cylindrical geome¬ 
tries, giving novel equations for the defect case, and conhrming also the previously 
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conjectured transmission factors. 


4.1 Lee-Yang Lattice Model and Transfer Matrices 

The Lee-Yang lattice model is defined on a square lattice with spins or heights a = 
1,2, 3,4 restricted so that nearest neighbor heights differ by ±1. The spins thus live 
on the A 4 Dynkin diagram. It is helpful to regard the Lee-Yang model as a special 
case of the Al Forrester-Baxter |49) model. 


4.1.1 Lee-Yang lattice model as the A 4 BF model 


The Al Forrester-Baxter ra models with spins a = 1,...,L are dehned by the 
Boltzmann weights 


W 


W 


W 


/ 

a ± 1 

a 

\ 




u 

V 

a 

a =F 1 


/ 

a 

a ± 1 

\ 




u 

V 

a =F 1 

a 


/ 

a 

a ± 1 

\ 




u 

V 

a ± 1 

a 

/ 


s(A — u) 

s(A) 


9aTi f s((a±l)A) \ s(u) 
ga±i V ■s(aA) J s(A) 


s(aA ± u) 

s(aA) 


(4.1) 


(4.2) 


(4.3) 


Here we differentiate between the massive and the massless cases as follows: 

• s(u) = sinu , the trigonometric sine function for the massless model 

• s(u) = 'di(M, q) the standard elliptic theta function for the massive model where 


i9i(u, q) = 2 g^'^^sinM ]^(1 — 2 ^^"" cos 2 m -|- q^AA ~ (4-4) 

n=l 

where u is the spectral parameter and the elliptic nome g is a temperature-like variable 
corresponding to the (^ 1^3 integrable bulk perturbation. 
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In our massive calculations for the Lee-Yang model we have found that C = 
where 

^=II}L= lim Nf 

4v3 N^co,q^0 

or more precisely 


L = lim Ya, 

A^^OO, a^O 


, 4v^d^ 

m = hm - 

—^0, Q-—^0 Ci 


where a is the lattice spacing, m is a mass, L is the continuum length scale, and t = 
is the deviation from critical temperature variable, and v = ^. 

The crossing parameter A is 

^ ^ {p' -p)t^ 

p/ 

where p < p'. Integrability derives from the fact that these local face weights satisfy 
the Yang-Baxter equation. 

Here we will only consider the Lee-Yang model with 



p' = 5, p = 2, X=- 

The gauge factors pa are arbitrary but we will take Pa = (—1)“'^^, so that = —1. 
With this choice the face weights are only symmetric under reflections about one 
of the two diagonals. Since the Lee-Yang model is non-unitary [p' ^ p + 1), some 
Boltzmann weights are negative. 

The critical theory correspond to r = 0 and then the function s{u) degenerates simply 
to s{u) = sin(M). 


4.1.2 Transfer matrices and functional relations 

The transfer matrices Ti{u) can be built up from the local face weights. As the 
local face weights satisfy the Yang-Baxter equations the transfer matrices will form 
commuting families [Tj{u),Tj/{u')] = 0 from which integrable Hamiltonians can be 
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derived. Due to the underlying symmetry of the model the transfer matrices will 
satisfy the functional relations 

Ti(n)Ti(n + A) =I + T 2 (n + 3A) (4.6) 

where in our case T 2 {u) is spectrally equivalent to Ti{u). Let us see how this functional 
relations are realized in the various circumstances. 


4.1.2.1 Periodic boundary condition 


The transfer matrix with periodical boundary condition can be dehned on a lattice 
of 2N sites from the local face weights as follows: 


Tj{ut = W, 


^ bi 62 

\ 


f 62 bs 

\ 

( 


u 

W, 


u 

...IT, 

^ ai 02 

) 


y 0-2 03 

) 

V 


&2V-1 

a2Ar-l 


&27V 

\ 


f b2N bi 


u 

IT, 


a2N 

) 


1 a2N Oi 


and W = Wi- This transfer matrix is a hrst in a series of transfer matrices obtained 
by fused weights. The transfer matrix T 2 {u), for instance, is dehned from the fused 
face weight W 2 : 


[ d c 

\ 

u 


^ e / 

\ 

u 

W 

d c 

\ 

n -|- A 

\ a b 

) 

^ ^ e 

^ a h 



V " ^ 

/ 


being non-vanishing only when \a — d\ = {0, 2} and a -|- d = {4, 6 }. It is not hard to 
see that the dehnition is good, i.e. the fused weights does not depend on /. These 
transfer matrices form a simple fusion hierarchy: 


Ti(n)Ti(n + A) 


s{u + A)s(A — u) 


N 

1 + 



N 

T,(u) 


Here N is taken to be even. The Lee-Yang theory is the simplest theory as in this case 
the fused weights W 2 are trivially related to Wi modulo some n-independent gauge 
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factors: 



(, 

\ 


( , 

\ 


d c 


d c 

W 2 


u 

= W 


M + 3A 


^ a b 

) 


^ 5 — a 5 — 6 

/ 


This gauge transformation can be dropped at the level of the transfer matrix, thus by 
introducing the height reversal matrix Y we conclude that Ti{u) = YT 2 {u) = T 2 {u)Y. 
If we renormalize the transfer matrix as 


T{u) 


( s(\Hu + 2\) 

+ A)s(u + 3A) / ^ 


then T{u) will satisfy the functional relation 


T{u)T{u + \) = l + T{u + 3\)Y 


As the height reversal matrix Y commutes with T we can diagonalize it in the same 
basis. Since Y"^ = 1 the eigenvalues are Y = ±1. Restricting the analysis to the 
Y = +1 eigenspace the eigenvalues of the transfer matrix t{u) will satisfy the relation 


t{u)t{u + A) = 1 + t{u + 3A) 

Clearly for any eigenvalue —t{u) is an eigenvalue, too. 
The transfer matrix also satisfies the crossing relation; 


T(A -u)= T{u) 


and periodicity 


T{u + tt) = T{u) 
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4.1.2.2 Periodic boundary condition with a seam 

The transfer matrix with periodical boundary condition with a seam of parameter ^ 
can be dehned on a lattice of iV — 1 (with even N) sites from the local face weights 
as follows: 


Tjiu)l = W, 



u + nwj 



u\ ...Wj 


b2N-l 

b2N 

a2N-i 

a2N 

\ N-1 

) 

IS a 


u I Wj 



factor which ensures that the transfer matrices satishes 


Ti{u)T\{u + A) — 1 + T 2 {u + 3A) 

where T 2 {u) is dehned from the fused face weights and fused seam. Just as in the 
periodical case we have Ti{u) = T2{u). We note that the ^ t 0 limit reproduces the 
periodic result. 

4.1.2.3 Boundary case: donble row transfer matrices 

To ensure integrability of Lee-Yang lattice model in the presence of a boundary m we 
need commuting double row transfer matrices and triangle boundary weights which 
satisfy the boundary Yang-Baxter equations. The integrable boundary conditions are 
labeled by the Kac labels (r, s) = (l,s) with s = 1,2. For s = 1, the non-zero left 
and right triangle weights are given by 


Kl 

c 

2 

u 

s(2A) 

Kr 

( 1 

2 

\ 

u 


Y 

) 

s(A) 


/ 


1 


(4.7) 
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For s = 2, the non-zero right boundary weights can be obtained by placing a seam in 
front of the s = 1 boundary 


Kr 


1 


^ a 2 


(2 1 


1 


2 

«,{ ^ W 


u + nw 


+ C Kr 

2 

u 


) 

\ 2 1 

J 

^ 0 2 

J 


J 


(4.8) 


s{u ^ (2 — a)X)s{u — ^ + (2 — a)X) 

s{X)s{2X) 


0 = 1,3 


The parameter ^ is arbitrary and can be complex. However, to obtain conformal 
boundary conditions at the isotropic point u = A/2, we choose 3fJe(^) = 0. Integra- 
bility in the presence of these boundaries derives from the fact that these boundary 
weights satisfy the left and right boundary Yang-Baxter equations respectively. 

The face and triangle boundary weights are used to construct [731 ^ family of com¬ 
muting double row transfer matrices D{u). For a lattice of width N, the entries of 
D{u) are given diagrammatically by 


i2(u)a.b = 5] 


Kl 


Co Cl 
1 ai 

X-u\w 


u W 


Cl C2 


ai 02 


u\ ...W 


CJV-I s 


a-N-i s 


u\ Kji\ cn 


U X 


1 hi 


Co Cl 


X-u\W 


hi 62 


Cl C2 


A-m ...W 


5jv-i s 


cn-1 s 


X — u 


This transfer matrix is positive dehnite and satishes crossing symmetry D{u) = 
D{X — u). Although D{u) is not symmetric or normal, this one-parameter family of 
transfer matrices can be diagonalized because L){u) = GD{u) = D{u)'^ is symmetric 
where the diagonal gauge matrix G is given by 


N-l 

Gab G{cij^ 5{cij^hj^ (4.9) 

1=1 

with 


G(a, 6 ) 


s(A)/s(2A), b 


1,4 


1 , 


otherwise 
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It is convenient to define the normalized transfer matrix 


T(m) = Ss{u) 


s‘^{2u — A) 


s{2u + X)s{2u — 3A) Ls(m + X)s{u + 3A)- 


s{X)s{u + 2A) 


N 


D{u) 


(4.10) 


with 


Ss{u) = 


1 , 


s = 1 


s^{X)s{u + ^ + 2A)s(m — ^ + 2A) 


(4.11) 


s{u + ^ + X)s{u — ^ + X)s{u + ^ + 3A)s(m — ^ + 3A) ’ 


s = 2 


It can then be shown pa that the normalized transfer matrix satishes the universal 
TBA functional equation 


T(m)T(m + A) = 1 + T(m + 3A) (4.12) 

independent of the boundary conditions. Apart from the change in the value of 
the crossing parameter A, it is also the same functional equation that holds for the 
tricritical hard square and hard hexagon models [68]. However, this change in the 
crossing parameter drastically changes the relevant analytical properties. Since the 
transfer matrices commute this functional equation also holds for each eigenvalue t{u) 
olTiu). 

4.2 Classification of Excited States 

Here we start the classihcation of states in the critical case when s{u) = sin(M). We 
will make correspondence with the conformal Lee-Yang model thus we recall first the 
description of their Virasoro modules. The Virasoro algebra contains two relevant 
modules for c = — ^. 

The identity module, is built over the vacuum |0) 

^ |0) ; Urn > 1 ; ni> Ui+i P 1 
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Interestingly this basis is linearly independent (no singular vectors). The representa¬ 
tion has the reduced character 


Xo(q) = ^dini(l^„°)g’^ = ^ 

n=l n=l 


q 






n 

n=l 


1 

(l-g5n-3)(l_^5n-2) 


where n can be considered as the energy (Lq eigenvalue) of the given state. The 
sum and the product form is related by the Andrews-Gordon identity, which is the 
generalization of the Rogers-Ramanujan identities. 

The other appearing module is built over the highest weight state \h) where h = — |. 
The module is generated by the linearly independent modes 


L—ni ■ ■ ■ ^ 0 j ^ ^i+1 T 1 


and has the reduced character: 


Xi{q) = 5^dim(lG)g’^ = ^ 

n=l n=l 


( 1(1 - 5 ”) 


n 

n=l 


1 

(1 - - g5n-l) 


4.2.1 {m,n) systems, zeros, paths and characters 

The Hilbert space of the lattice model consist of paths. By diagonalizing the various 
transfer matrices we can characterize a given eigenvector by the distribution of the 
zeros of the transfer matrix on the complex plane. We will make correspondence 
between the paths, the distribution of zeros and the Virasoro descendants in the 
three cases. We start with the simplest boundary case as the Hilbert space contains 
one single Virasoro module only. Then we turn to the analysis of the periodic case 
with and without the seam, where tensor products of Virasoro modules appear. We 
also analyze the flows between the Hilbert spaces induced by ^ going from 0 to cx). 
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4.2.1.1 Boundary case 

(m, n) system, characterization of the eigenvectors by zeros of the transfer 
matrix 

Let us consider the sectors with boundary conditions {r,s) = (1,1), (1,2) which we 
often label simply by s = 1, 2. The excitation energies are given by the scaling limit of 
the eigenvalues of the double-row transfer matrix D{u), or equivalently the normalized 
transfer matrix T{u), where u is the spectral parameter. The single relevant analytical 
strip in the complex u-plane is the full periodicity strip 

- I < Lte{u) < ^ (4.13) 

5 5 

Additionally, the transfer matrix is symmetric for the real line thus it is enough to 
analyze the analytical structure on the upper half plane. The excitations are classihed 
by the string content in this analytical strip. There are four kinds of strings which we 
call “1-strings", “short 2-strings", “long 2-strings" and “real 2-strings". See the next 
figure for two typical configurations in the two sectors. 


■ 

1-string 


• • 


• • 

• • 

J k 

short 2-string 

) k 

• • 

) ( 

) long 2-string ( 

) C 

—1-1-1-h-•- 

_ real 2-string _ 

1-1-d-h-•— 


5555 5555 


Figure 4.1; Typical conhguration of zeros of the transfer matrix eigenvalue. The 
boundary condition (1,1) is on the left, while (1,2) is on the right. 

The 1-string Uj = 37r/10 -|- ivj lies in the middle of the analytical strip and has real 
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part Stt/IO and exist in the (1,1) sector only. The two zeros of a short 2-string 
Uj = 7r/5 -|- iwj, 27r/5 -|- iwj have common imaginary parts and real parts 7r/5, 27r/5 
respectively. The two zeros of a long 2-string Uj = —7r/5 -f- iWj, dvr/S -f- iWj have 
common imaginary parts and real parts —7^/5, dvr/S respectively so that these zeros 
sit at the edge of the analytical strip. Lastly, a real 2-string consists of a pair of zeros 
Uj = 7r/10,7r/2 on the real axis. The string contents satisfy the {m,n) system 

2m + n + 3 — s = N^m + n = N — 3 + s, s = l,2 (4-14) 

m = {number of short 2-strings}, n = {number of long 2-strings} (4.15) 

There is always a real 2-string on the real axis and, in the (r, s) = (1,1) sector, a 
single 1-string furthest from the real axis. Each “short 2-string" contributes two zeros 
and, by periodicity, each “long 2-string contributes one zero. The 1-string contributes 
one zero and so does the real 2-string since it is shared between the upper and lower 
half planes. Consequently, the (m, n) system expresses the conservation of the total 
number of zeros in a periodicity strip. The roles of m and n are interchanged under 
duality. For the leading excitations m is hnite but n ^ N as N ^ oo. 

As explained in m, an excitation with string content (m, n) is uniquely labeled by a 
set of quantum numbers 

J = (/i,/2,...,Jj (4.16) 

where the integers Ij > 0 give the number of long 2-strings whose imaginary parts Wj 
are greater than that of the given short 2-string Vj. The short 2-strings Vj and long 
2-strings Wj labeled by j = 1 are closest to the real axis. The quantum numbers Ij 
satisfy 

n>Ii>l2>--->Im>0- (4.17) 

For given string content {m,n), the lowest excitation occurs when all of the short 
2-strings are further out from the real axis than all of the long 2-strings. In this case 
all of the quantum numbers vanish Ij = 0. Bringing the location of a short 2-string 
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closer to the real axis by interchanging the location of the short 2-string with a long 
2-string increments its quantum number by one unit and increases the energy. 


Finitized characters 

For N = s — 1 (mod 2), the (fermionic) hnitized characters are 


(r,s) = (1,1) : 

(r,s) = (l,2): 


(7V-2)/2 


xZ= E 


m=0 


(7V-l)/2 


N — 2 — m 


m 


Xo{q) 


-I Q 


m=0 


iV — 1 — m 
m 


Xi{q) 


(4.18) 

(4.19) 


where 


N 

m 


n 


l-q 


N+l-i 




n 


=1 ^ i=i ^ 


for iV —)■ cx) 


These hnitized characters can also be written in the form 


CZ = 


E«'" = E« 




.o-j + i) 


(4.20) 


where the sum is over all one-dimensional RSOS paths a = {aQ,( 7 i,... ,a]\f} on A 4 
with (To = s and (Ttv = 1- The energy function H is 


H(^(jj—i, (jj, (jj_|_i) 


I 1, (^j, (^j+i) = (2,1, 2) or (3,4, 3) 

0, otherwise 


(4.21) 


Notice that this local energy function differs from the one introduced by Forrester 


and Baxter. 
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Bijection between RSOS paths, strings and Virasoro modes 

There is in fact a bijection [50] between the one-dimensional RSOS paths a that label 
the eigenstates (eigenvalnes), the allowed patterns of strings in the periodicity strip 
and the state described in terms of the Virasoro modes. A triple (cxj-i, Uj, Uj+i) = 
(2,1,2) or (3,4,3) corresponds to a short 2-string (particle) at position j and an 
insertion of a Virasoro mode whereas a pair segment (cTj,cTj+i) = (2,3) or (3,2) 
corresponds to a long 2-string (dual particle) at position j + ^- In addition, in the 
sector (r,s) = (1,1), there is a 1-string at j = 0 corresponding to the initial height 
s = 1 at j = 0. This bijection is illustrated in Figure [T^ Notice that only the relative 
positions of the long and short 2-strings is important. If s = 1 the hrst and last (3 —s) 
segments are inactive whereas, if s = 2, only the last (3 — s) segments are inactive. 
We see that the geometric constraint 

2m + n + 3 — s = N (4.22) 


agrees with the (m, n) system. 

Flow between boundary conditions 

We are in the position now to describe the boundary flows induced by ^ between the 
boundary conditions (1,1) and (1, 2). This flow is realized when ^ goes from 0 to cxd 
and we can describe it at the three different languages we already introduced. 

In terms of the zeros and (n, m) system the flow is very simple: the 1-string which 
exist only for the (1,1) boundary condition start to move to inhnity in the imaginary 
direction as indicated on the hgure. There is no change in the 2-strings. 
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Figure 4.2: Bijection of one-dimensional configurational paths a, strings in the peri¬ 
odicity strip for the (r, s) = (1,1) sector and the Virasoro modes over |0) is shown on 
the left part. The right part contains the analogous hgures for the case (r, s) = (1, 2). 
The paths a (rotated 90° clockwise) are shown on the left and the analytical strip 
containing a 1-string (solid square), short 2-strings (solid circles) and 2-strings (open 
circles) in the upper-half complex u-plane are shown on the right. For each short two 
strings we associate a Virasoro mode where they are located, L _2 and L _5 in the (1,1) 
and L_i , L _4 in the (1,2) cases. The string (particle) content is m = n = 2 with 
iV = 8 in the (1,1) and with iV = 7 in the (1, 2) cases. The energies of these excited 
states are E = 2 + 5 = 7 (i? = 1-1-4 = 5) since there are two particles (short 2-strings) 
at respective positions j = 2,5 (1,4), respectively. The corresponding states in the 
Virasoro modules are L_5L_2|0) and L_ 4 L_i|/i). 
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5555 5555 5 55 5 

Figure 4.3; The boundary flow in terms of the zeros. On the left = 0 while on the 
right ^ = oo. In the ^ oo limit, shown in the middle, nothing else happens than the 
imaginary location of the 1-string moves to infinity and disappears. The distribution 
of the 2-strings are not affected at all. 

The flow in terms of the paths is also simple we merely have to remove the first raw 
of the pathspace. 

Most enlightening is the flow in terms of the Virasoro modes. First of all the highest 
weight state |0) flows to |<h), and in the module the rule is very simple we have to 
increase the index of every Virasoro mode by one —>■ L_„+i: 

L—ni L—n2 ■ • Hfc |0>^L —n\+lL—n2+l---L—nk+l\h) 

This very simple flow is summarized for the first few excited states in the following 


table 
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Level 

State in the (1,1) module 

State in the (1, 2) module 

Level 

h.w. state 

|0) 

|<I>) 

h.w. state 

2 

L-2i0} 

L-im 

1 

3 

L-3i0} 

L-2m 

2 

4 

L-M 

L-3i^} 

3 

5 

L-5i0} 

L-4m 

4 

6 

-h-2-h-4 0) 

L-iL.3m 

4 

6 

L-eiO) 


5 

7 

-h-2-h-5 0) 

L-iL_4j^} 

5 

7 

L-rjO) 

L_6|<h) 

6 

8 

-h-2-h-6 0) 


6 

8 

-h-3-h-5 0) 

-h-2-h-4 4*) 

6 

8 

i^-8|0) 

L-rm 

7 

9 

-h-2-h-7 0) 

L-iL_ej^} 

7 

9 

-h-3-h-6 0) 

L-2L.5m 

7 

9 

L-gjO) 

L_8|4>) 

8 






As expected and shown in the table above, the character will flow from Xi^i = l + g^ + 
g^ + g'^ + g® + 2g® + 2g’^ + 3g® + 3g® + ... to Xi ^2 = l + g + g^ + g^ + 2g'^ + 2g^ + 3g® + 3g'^ + .... 

The level by level flow agrees with the TCSA resnlt of jlHI 

4.2.1.2 Periodic case 

The best way to describe the periodic case is based on the previously introduced 
boundary identification. In analyzing the zeros of the transfer matrix we can distin¬ 
guish two different appearance of zeros as shown on the figure 4.4 
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JL 23 
5 5 


23 

5 


of 


Figure 4.4; The two appearing configurations of zeros of the transfer matrix eigenvalue 
in the periodical case. 


The first we can recognize is that we have similar short and long 2-strings and 1-strings 
as in the boundary case. What is different, however, is that the zeros on the lower 
half plane are not necessarily related to those on the upper half plane except for the 
1-string. So if we have a 1-string on the lower half we always have one the upper half, 
too. In classifying the states we can use the already developed classification for the 
boundary case, taking into account that the lower and upper halves are independent. 

For the (m, n) structure, we have to differentiate between the structures on the 2 
sides of the real axis. Now our lattice is 2N dimensional with N zeroes on each 
side. We define an (m, n; fh, n) system. On each side the eigenvalues with zeros at || 
correspond to the s = 1 and have N — 2 zeros, while the other eigenvalues he in the 
s = 2 sector with N zeros. Formulating this, we get that 2m + n — 2s-\-4 = N or 
equivalently 2m + n = N + 2s — 4, where m is the number of the short 2-strings and 
n is the number of long 2-strings. Similarly 2fh + n = N + 2s — 4. 

Studying the zero structures, we find that they resemble the Hilbert spaces of V 0 Vi 
for the structures with zeroes with the 1-string and V 2 0 V 2 for those without the 
1-string. In fact, the real axis separates the Vi and V 2 parts of the tensored states 
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Figure 4.5; The identification between the zeros of the transfer matrix, paths and 
Virasoro modes in the periodical case. 


from the Vi and V 2 respectively. We define the vacuum |0) as the state with |0,0) and 
the (1,2) ground state |<F) as the |^, state. 

The state with one short string furthest from the real axis up below the 1 string 
corresponds to T_ 2 | 0 ) and moving the short string downwards through the long strings 
increases the level by 1 for each permutation, thus creating the L_„|0) states. The 
mirror image on the zeroes below the real axis corresponds to the L_n\U)- A similar 
description also applies for L_„j...L_„^|0). 

For the (1,2) states, the lowest excitation appears with a short string on the top of 
all long strings, with no 1-string above. This is L_i|<h), and every time we lower the 
short string below a long string we obtain one extra unit of energy, hence we have all 
the L-n\^) and similarly for the mirror image of L and for combinations of those. 
Summarizing, in this model we find out that the first few states that we obtain from 
the classification of the zeros of the eigenvalues are: 
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Level 

Fi (g)Fi 

V 2 <S> V 2 

0 

|0) 

m 

1 



2 

L_2|0), Z_2|0) 

l_2|$), Z-2|$), L-iL-i|1>) 

3 

l_3|0) , Z_3|0) 

l_3|$), L-3|$), L-2L-i|$), i-lL-2|$) 

4 

L_4|0) , L_4|0), L_2L_2\0) 

-^- 4 !^)) L_3L_i|<t>), L_2L_2\^), 



Z_3L_i|4>), Z_ 3 Z_i|$), and Z_ 4 |$) 


Table 2 lists the first few states that we can see from the zero eigenvalues of the 
periodic transfer matrix. In a lattice of 2N sites, one can read all the states up to 
level N completely. The same can be done from the RSOS paths in the Hilbert space. 

4.2.1.3 The case of a seam 

Introducing a seam we can analyze the two limiting cases similarly we analyzed in the 
boundary setting, namely ^ going from 0 to cx). Clearly for = 0 the seam disappears 
(identity seam) and we recover the results of the periodic boundary condition. In the 
^—>■00 limit we found the following identihcation between the strings the paths and 
the Virasoro modes 



Figure 4.6; The identihcation between the zeros of the transfer matrix, paths and 
Virasoro modes in the case of a seam at the hrst site. The seam shadows the periodicity 
of the paths. 
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The flows are very simple in terms of the zeros. We have the following simple mech¬ 
anism. The flows can be explained in three mechanisms: 

A. If the outermost string is a 1-string, it flows towards infinity with increasing b. 
(Plot on the left) 

B. If the outermost string is a short 2-string, one of the zeroes flows to infinity and 
the other goes to (Plot on the right) 

C. If the outermost string is a long 2-string, it flows towards infinity. (Plot in the 
middle) 

In terms of the states this is summarized as follows: 

1. Due to type A flows: 

with |0) —)■ 10) 

2. Due to type B flows: 

3. Due to type C flows: 

Now using the first few states of the trivial defect case from Table 2, we will deduce 
their corresponding states using this mechanism, and it is indeed what we can observe 
from the flows of the zero eigenvalues as was shown for sample states above. 
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Level 

Trivial Defect 

Non-trivial Defect 

Level 

h.w. state 

|0) 

10) 

h.w. state 

h.w. state 

1$) 

1$) 

h.w. state 

1 


10) 

h.w. state 

1 

T-i|4>) 

T-i|4>) 

1 

2 

T- 2 IO) 

Z- 1 I 0 ) 

1 

2 

Z_2|4>) 

Z-i|4>) 

1 

2 

T-iL_i|$) 

L-M 

1 

2 

L_2|4>) 

L_2|4>) 

2 

2 

L-2\0) 

T_2 0) 

2 

3 

L-3\0) 

Z- 2 I 0 ) 

2 

3 

L_3|4>) 

L_2|<h) 

2 

3 

T-2T_i|4>) 

Z-iL_i|$) 

2 

3 

T-iL_2|4>) 

L-M 

2 

3 

L_3|4>) 

L_3|4>) 

3 

3 

L-3\0) 

^-3|0) 

3 

4 

Z_iZ_3|4>) 

Z- 2 I 0 ) 

2 


Table 3 shows the exact flow from each state in the trivial defect Hilbert Space to its 
corresponding state in the non-trivial one up to the second order descendent level in 
the defect Hilbert space. 

It can be seen that the flow occurs from Hi 0 Id +® V 2 to Id ® Id + Id ® Id + V 2 ® Id 
as we expect from defect conformal field theory. Had we taken the limit b —)■ — 00 , 
we would have got the similar total outcome, but with 0 and 0 exchanged, and 
being the operator augmented to instead of L as happens here. 
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Figure 4.7: The identification between the zeros of the transfer matrix, paths and 
Virasoro modes in the case of a seam at the first site. The seam shadows the pe¬ 
riodicity of the paths. The two figures on the left correspond to type A flow, the 
outermost 1-string flows to infinity, and the state flows from L_ 3 L_^L_ 2 L_^\ 0 ) —)■ 
L_3L_5L_iL_ 4|0). The middle two figures are of type C, the long 2-string goes 
away to infinity and the state flows from L_ 4 L_ 6 T_ 2 T- 5 |‘F) —?• T_ 4 L_ 6 T-iT_ 4 |<h). 
Finally the two figures on the right show how a short 2-string has one of its ze¬ 
roes flow to infinity and the other to the middle n = |^, and the state flows from 

T_4L_6T_iZ_5|$) —)■ L_4L_6T_4|0) 


4.3 TBA Equations 

In this section we solve the TBA equations of the Lee-Yang model for the periodic 
boundary conditions with and without a seam, and for the boundary model in both 
the critical and the massive cases. We derive those equations on the lattice, and after 
scaling we confirm the results of the continuum limit equations which were derived 
in inn EH inni. Our approach is systematic since we know the analytic structure and 
the zero eigenvalue locations in the analytic strip, which allows us to solve the TBAs 
taking into account this structure. 
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4.3.1 Critical/Massless TBAs: 

4.3.1.1 Periodic boundary conditions 

The transfer matrix satisfies the functional relation 


T{u)T{u + A) 


sin(A + u) sin(A — u) 
sin^ A 


N 


+ 


Sinn 

sin A 


N 

T{u + 3A) 


We normalize the transfer matrix as: 


t{u) 


( s{\)s{u + 2\) y 

\s(n + \)s{u + 3A )) ^ 


and we get that 


t{u)t{u + A) = 1 + t{u + 3A) 


Using the periodicity t{u) = t{u + tt) we rewrite it as t{u)t{u — ^) = 1 + t{u — 
and after shifting u we have 


, TT, , TT, ^ , , 

t{u + -)t(n --) = ! + t{u) 

5 5 

We decompose t into two components / and 1: 

t{u) = f^{u)l{u) 

where corresponds to the bulk free energy (order N term) and I corresponds to 
the hnite size corrections (order A term) and N is even. We want to kill the iVth 
order zeros at — ^ and poles at ^ in /(n) which satishes 


oi|ii 
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The solution compatible with the analytical structure is 


/(w) 


sin ^ + sin I 
sin - sin I 


,5m tt, ,5m 

= - tan(— + -) tan(— + 2-, 
6 6 6 6 


which is basically the shifted S-matrix. /(m) satishes /(m)/(m + A) = 1. Introducing 
the variable m = ^ we write the functional equation as 


t{x 



)t{x + i^) 


1 + t{x) 


In this variable 


., , Ax 57r, Ax Ttt, 

f{x) = - tan(y + —) tan(— + —) 


(4.23) 


Vacuum state: 

The ground state has no zeroes inside the analytical strip, and since it is analytic in 
the strip, we can use the functional relation to write 


t{x — i^)t{x-L i^) l{x — i^)l{x + i^) ^ ^ 


t{x) 


l{x) 


= l + t (x) 


For the ground state, inside the physical strip x E y]) both sides of the equation 
are ANZ, so we can take the log and solve in Fourier space, then 


\ogl{x) =-(f-k\og{l + t \x)) ; Lp = - - ^^ 

1 — 63—e 3 


where 


and 


zn / 


(p-kf= / dy^{x-y)f{y) 


—oo 


(4.24) 
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Now we restore t: 


\ogt{x) = N\ogf{x) — ip'klogi^l + t ^(x)) (4.25) 

This is the ground-state TBA on the lattice. In the thermodynamic limit all inter¬ 
esting things happen around two domains; in the u variable either on the upper half 
plane around |ilogiV or on the lower half plane around — logiV. For this reason 
in the x variable we focus on the behavior around ± log N. Let us center the new 
functions around x ± logiV as = lim t{x ± logiV). Taking the continuum limit 

N^oo 

{N —)■ oo) on the source term we get 

TT 1 TT 

N lim log/(a; -|- logiV) = N lim log(l -|- 4sin —e^^— -!-...)= 4sin 
N^oo N^oo 3 N 3 

which leads to the massless ground-state TBA equations 

€^{x) = 4 sin — ific log(l 4- (4.26) 

o 

where 4 sin | can be easily absorbed by shifting x. 

Excited states: 

From our numerics, we know that in the excited states there are zero eigenvalues in 
the analycity strip. We classihed them before as 1-strings and short 2-strings. The 
1-strings occur at 

Stt 

u = — + za 

and additionally the short strings at 



u = 
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For finite energy states in the continuum [N —)■ cx)) limit, they go to inhnity as 
a = |(±logiV + a^) and I3j = |(±logiV + on the upper/lower half plane, 
respectively as it can be analyzed numerically. 

In the X = variable they are located at 


Xg = ^ = log iV + q;=‘= 
xf = ±i| + log iV + j3f 

As we would like to take logarithm we need functions free of zeros and poles on the 
physical strip. The function which can eliminate the single zero is 


(To(a;) = tanh 



while the one which eliminates the two zeros at ±f is 

D 


ai{x) = f{x) ^ 


cosh a; + sin | 
cosh X — sin | 


These functions satisfy 


aoix - f)ao{x + ^) = 1 
ai{x - f)(Ti(x + f) = cri{x) 


(4.27) 


(4.28) 


(4.29) 


With these functions we parametrize the normalized transfer matrix eigenvalue as 

M 


t{x) = f{x)^ JJ ao{x - x^) JJ cTi(a; - xf)l{x) (4.30) 

± i=i 

which ensures that l{x) is ANZ in the physical strip. The functional equation then 
takes the form 


JJ(To(a; - 

± 



t{x — i^)t{x + f |) 
t{x) 


l{x — i^)l{x + f |) 
l{x) 


JJao(a;-a;^)(l + t ^(x)) 

(4.31) 
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Clearly both sides are ANZ in the interior of the physical strip; the combination 
n± o'oix — XQ)t~^{x) is regular at a; = Xq. Taking then logarithm and going to Fourier 
space we hnd; 

log/(a;) = -(/?^log(JJao(a; + f^(a;))) 

± 

Restoring t we obtain the excited-state massless TEA equation on the lattice: 

logt(a;) = A^log/(a;) + log(ao(a; - log(cri(a; - xf)) 

-ip-k\og(Jl^ao{x - Xo){l + t-\x))) 

The parameters of the excited state Xi = {x'^,x^} are determined self-consistently 
from the fact that it is a zero of the transfer matrix; 

^(a^)L=xi±f = -1 

In the scaling limit we can focus on the two scaling domains at ± log N by introducing 
qUO = lim (To(a; ± log A^ — Xg -I-log A^) 

N^oo 

and in this limit, it satishes the excited-state massless TEA equation; 


e^(a;) = 4sin + ^logai(a; -- 9 ?^ log(ao(a; - (5=^)-h e (4.33) 


The parameters satisfy the following equations 




4.3.1.2 Seam 

The transfer matrix of the model with periodical boundary conditions with a seam 
of parameter ^ can be dehned on a lattice of N sites (A^ is even) from the local face 
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weights and satisfies the relation; 


TAu)TAu+X) 


s{u + A + ~ u — 0 f s{u + A)s(A — u) 

s2(^) ly s2(X) 


N-1 


1 + 


s{u + 0 
s{X) 



N-1 

TAu+3X) 


Defining the normalization 


_ s{\)s{u + 2A + ,^) / g(A)g(-u + 2A) ^ 

s{u + A + ^)s{u + 3A + .^) \s{u + X)s{u + 3A)/ ^ 

We obtain the functional relations; 


t{u)t{u + A) = 1 + t{u + 3A) 

Using the same relation assumed for the periodic boundary conditions, we would like 
to kill the (iV — l)th order zeros at —^ and poles at and the order one zeroes 

at—| —and poles f — ■C) Here ^ = ib, is a pure imaginary parameter. 

We use the same 

, ,5m tt, ,5m tt, 

/(m) = - tan(y + -) tan(y + 2-) 

and write the transfer matrix in the form 

t{u) = f{u)^-^f{u + Ol{u) 

We introduce the variable m = ^ and solve functional equation 

t(x - i^)t{x + i^) = 1 + t{x) 
o o 


Vacuum state 


Since this state has no zeros in the physical strip we get that; 


logt(a;) = (iV-l)log/(a;)+ log/(a;--(^^log(l + t ^(x)) (4.34) 
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This is the ground-state TBA on the lattice with a seam. 

In the thermodynamic limit u scales as ±z| logiV (upper or lower half plane), hence 
in the x variable we focus on the behavior around ± log N. We can also play with the 
parameter If we do not scale it in the thermodynamic limit it simply disappears 
from the equations. If we scale with =pi|logA^ it will appear in the equation for 
only, respectively. Let us focus on the plus sign: so ^ = —^|(^ + logA^). Now we 
center the new functions around x ± log as = lim t{x ± logiV). Taking the 

N^oo 

continuum limit (TV —)■ oo) on the source term we get the massless ground-state TBA 
equations in the presence of a seam; 

e^(x) = 4 sin + \ogf{x — 0 — v^T^log(l -|- (4.35) 

o 

This agrees with the scaled bulk TBA equation as 

/(x-f) =T_(y+a;) ; 6 = 3+^ 

2 71 

where T_{x) = [-1 + b][l + b] with [s] = 

T_ is the transmission matrix in the continuum theory with a defect for b = 3 + ^ 
as we expect from |37] . 


Excited states 


For large where the flow of states has already occurred to the (1,2) module, we 
might have 1-strings at 


u = 


Stt 

10 


la 


and short 2-strings at 



u = 
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For finite energy states in the continuum [N —)■ cx)) limit they go to inhnity as 
a = |(±logiV + Q;^) and fij = |(±logiV + /3j*"), but here a occurs either in the upper 
or in the lower half plane, and (if need not be symmetric with respect to the real 
axis. 

In the a; = — 1^) variable they are located at 

= ^ = logiV + d± 
xf = ±i| + log iV + (3f 

Again, it is important to note that the ± sign here on x is to indicate whether the 
zero is occurring in the upper or the lower half plane, and doesn’t indicate symmetry 
with respect to the real axis. Similar to the periodic case, to take logarithm we 
need functions free of zeros and poles on the physical strip. The function which can 
eliminate the single zero is cro(x) = tanh (|a;) while the one which eliminates the 
two zeros at ±f is ai(x) = f(x)-^ = • 

With these functions we parametrize the normalized transfer matrix eigenvalue as 

M 

t{x) = f{xYaQ{x-x^)W_ai{x-xf)l{x) (4.36) 

i=i 

Here, there is no product over ao as it only occurs once, either in the upper or in 
the lower half planes. With this parametrization, l{x) is ANZ in the physical strip. 
Following the same derivation of the excited periodic states, we hnd; 

log/(a;) = -(^^log(cTo(a; + f^(a;))) 

Restoring t we obtain 

logt(a;) = (iV - 1) log/(a;) + log/(a; - ^) + log(ao(a; - 
+ Ei log(c^i(a^ - xf)) -T* log(cTo(a; - Xq )(1 + f 


(4.37) 
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In the scaling limit we, we repeat same scalings as before and we get that = 

lim ao{x ± log— XQ)~^t{x + logA^) satishes the equation; 

N^oo 


e^{x) = 4sin + log/(x - 0 + J2j log<^ 1 ( 2 ^ “ Pf) 

—if * log(cTo(x — a^) + 


(4.38) 


where, as mentioned before f {x — ^) = T-{y + x) ; b = 3+^ and the parameters 
satisfy the following equations — Ci^)\x=xi±if = 


4.3.1.3 Boundary 

It is convenient to dehne the normalized transfer matrix for the boundary double row 
transfer matrix as 


t{u) = Ss{u) 


s^(2u — A) 


s{\)s{u + 2A) 
s{2u + A)s(2m — 3A) Ls(m + \)s{u + 3A) 


N 


D{u) 


and 


Ss{u) 


1 , s = 1 

s(A)^s(M+^+2A)s(ti—^+2A) _ ^ 

s{u+^-\-X)s{u+^+3\)s{u—^+3X)s{u—^+\) ^ 


where N is even for s = 1 and odd for s = 2. It can then be shown ca that the 
normalized transfer matrix satishes the universal TBA functional equation 


t{u)t{u + A) = 1 + t{u + 3A) 


(4.39) 


Using the periodicity t{u) = t{u + vr) we rewrite it as t{u)t{u — ^) = 1 + t{u — |) and 
after shifting u we have 

t{u + - ^) = 1 + t{u) 

5 5 

We write 

t{u) = f{uf^g{u)l{u) 
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where is the bulk free energy (order N term), g is the boundary energy (order 
1 term) and I corresponds to the hnite size corrections (order ^ term). Similarly as 
before we would like to kill the iVth order zeros at —^ and poles at ^ by f{u) 
which satishes /(« + |)/(n — |) = f{u) where the solution f{u) compatible with the 
analytical structure is 


/(«) 


sin ^ + sin I 
sin Y ~ sin | 


5m 71 5u 

tan(— + -)tan( — 

DO D 


+ 


2 % 

~6 


) 


which is basically the shifted S-matrix. 

Introducing the variable m = ^ we write the functional equation as 


t{x 



)t{x + i^) 


1 + t{x) 


(4.40) 


and in this variable, 


ix dvr ix Ttt 
f{x) = - tan(y + —) tan(— + —) 


(r,s) = (l,l) sector 

We also want to eliminate the zeroes of order 1 which originate from the normalization 
and from the analytical strip. 

For s = 1, Ss{u) = 1 , hence it has no contribution here. Hence the only contributions 
come from a double zero at m = | and poles at w = ^ = ^ and u = ^ = 

= ^. Due to I periodicity, it also has poles at | and Also from the analytic 
structure there are zeroes occurring at ^ and |. The factor that we need to insert 
and compatible with the analytical structure is: 

tan^(f - f) tan(f - p) tan(f - hp) 
tan(f - f)tan(f - f) 


91,lip) 
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In terms of the x variable we get: 


9iAx) = 


tan2(f)tan(f + f)tan(f - |) 
tan(f - tan(*f + 


where gi^i{x) satisfies the relation 


TT TT 

giAx - i-)gi,i{x + i-) = gi,i{x) (4.41) 

Here it is important to mention that gi^i{x) = + x)R^^{y ~ where Ri 

is the reflection matrix of the continnum boundary theory of the Lee-Yang model 
corresponding to the identity module, dehned by Ri{x) = (|)(|)(—2) and (s) = 
sinh(i-!^) • This will be important when we analyze the massive case. 

In addition, in this sector we always have one strings at 

Stt 

u = - ta 


and short strings at 


^ _ I I + 

For hnite energy states in the continuum [N —)■ oo) limit they go to inhnity as 
a = |(±logiV -|- and I3j = |(±logiV -|- I3f) on the upper/lower half plane, 
respectively. In the boundary case a~ = —a~^, and /3~ = as we can see from their 
symmetry with respect to the real axis. 

In the X = §i{u — ||) variable they are located at Xq = ^ = logiV -|- d^and at 
x^ = -|- logiV + j3f- As we would like to take logarithm we need functions free of 

zeros and poles on the physical strip. Those functions were dehned before as ctqIx) and 
(Ji{x). With these functions we parametrize the normalized transfer matrix eigenvalue 
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t{x) = f{xY^gi^i{x) JJ ao{x - Xq ) JJ ai(a; - xf)l{x) (4.42) 

± j=i 

which ensures that l{x) is ANZ in the physical strip. The functional equation then 
takes the form 


JJcro(a: - 

± 



t{x — i^)t{x + i|) 
t{x) 


l{x — i^)l{x + i|) 
l{x) 


+ f ^(x)) 


Clearly both sides are ANZ in the interior of the physical strip; the combination 
])([^(To(a: — x'^)t~^{x) is regular and non-zero at x = Xq . Taking the logarithm and 
going to Fourier space we hnd: 


log/(x) = -(p*log(JJcro(x -Xo)(l + t ^(x))) 

± 


where the convolution * was dehned in equation (4.24). Restoring t we obtain 


logt(x) = 2A^log/(x) +log^i,i(x) + ^^log(ao(x-x^)) + log(ai(x 

-p * iog(n± - xq )(i+ 


xf)) 

(4.43) 


The parameters are determined by t{x)\^^^.^in = —1. 

In the scaling limit we can focus on the two scaling domains at ± log N . We have 
calculated the scaling limits of /, cro(x — Xg )and cri(x — xf) before, with the new 
function (y'i,i(x) scaling to 1 in the massless case, hence disappearing from the scaling 
TEA. Using 

QpU) = lim cro(x ± log A^ — Xq )“^t(x-I-log A^) (4.44) 

V —>-00 


The parameters satisfy e'^*^®Vo(x — d^)| 


X=Xi±if 


- 1 . 
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And e{x) satisfies the equation 

e^{x) = 8 sin ^ log(Ti(a; — j3f) — 99 ^ log((Jo(a; — (4.45) 

3 

j 

This is the general massless TBA in the (1,1) sector. It is important to mention that 
for the groundstate of this sector, there are no short 2 -strings, hence Ui doesn’t appear 
and the ground state massless TBA of this sector on the lattice is: 


logt(a;) = 2iVlog/(a;) + log^i,i(a;) + ^^log(ao(a;-a;o)) 

-if -k log(n± )(1 + 

while the scaled ground-state massless TBA is 


(4.46) 


e^(a;) = 8 sin -I —99 ^ log((To(a; — d^)-|-e (4.47) 

o 


(r,s) = (l,2) sector 

Similar analysis follows in this sector, where we have to take into account the contri¬ 
bution of Ss{u). 

For s = 2 , S 2 {u) is proportional to f{u + ^)f{u — ^), therefore we should account for 
their respective zeroes and poles. 

The other contributions from the normalization and the analytic strip are the same 
as before, and they were included in gi^i{u). Then we need to dehne the new gi ^2 term 
as: 


9iM = 9i,iiu) X f{u + ^)f{u - 0 


In terms of the x variable we get: 


^1,2(3^) = 


tan 2 (f)tan(f + f)tan(f - f) 5iP 

■fix - ^)f{x +—, 


tan(y - f^) tan(Y -h A) 


(4.48) 
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where ( 71 , 2 ( 2 ^) satisfies the relation 

^ 1 , 2 ( 3 : - i^)9i,2{x + = 9i,2{x) 

Again we note that 9 i, 2 {x) = R'^^ +x) Ri^ ~ x) where Ri and R^ are the refiection 
matrices of the continuum boundary theory of the Lee-Yang model corresponding 
to the Id and (j) modules of the Virasoro algebra with highest weights 0 and ^ 
respectively. They are defined by 

fl* = (5)(i)(-2)(¥)(¥)(¥)(Y^) 

in the case where b = 3 — This will also be important in the analysis of the 
massive case. 

In addition, in this sector, in general we will have short 2-strings ( no 1-strings) located 
at; 

_ I I + 

In the X = ^{u — variable they are located at xf = ±i| -|- logiV -|- 
We parametrize the normalized transfer matrix eigenvalue as 

M 

t{x) = f{xf^9i^2{x) JJ cTi(a; - xf)l{x) (4.49) 

i=i 

which ensures that l{x) is ANZ in the physical strip. The functional equation then 
takes the form 

t{x - q)t{xq) i{x - q)i{x + i^) 

-fww-^ =-yww-1 + t X 

t[x) l[x) 

Both sides are ANZ in the interior of the physical strip. Taking then logarithm and 
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going to Fourier space we find; 

log/(a;) = —99*log(l + t~^{x)) 


Restoring t we obtain 

logt(a;) = 2iVlog/(a;) + log 5 fi, 2 (a^) + ^ log(o-i(a; - - 99^1og(l + t“^(a;)) (4.50) 

i,± 

The parameters are determined by = — 1 . 

^ 3 

In the scaling limit we again focus on the two scaling domains at ± log N. We have 
calculated the scaling limits of those functions before, with the new function ( 71 , 2 ( 2 ^) 
scaling to 1 in the massless case for hnite However, if we scale ^ —)■ —i|(.^ + logiV), 
and use = hmjv-s .00 t{x + logiV), it satishes the equation 

e^{x) = 8 sin + log 5 fi, 2 (a; - f) + '^\ogai{x - Pf) - 99 ^ 1 og(l + (4.51) 

o 

3 

This is the general TEA in the (1,2) sector. For the groundstate of this sector, there 
are no short 2-strings, hence ai doesn’t appear and the ground state massless TEA 
of this sector on the lattice is: 

logf(a;) = 2N\og f{x) + log ( 71 , 2 ( 0 ;) — v? ^ log(l -|- t~^{x)) (4.52) 

while the scaled ground-state massless TEA is 

e^(a;) = 8 sin -h log 5 fi, 2 (a; - i) - 99 ^ log(l -f- 
o 


(4.53) 
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4.3.2 Off-Critical/Massive TBAs: 

4.3.2.1 Periodic boundary conditions 

In the massive description we have to make the following replacements sinw —>■ 
i9i(u, q) and cosm —)■ 'd 2 (M, q). We solve the fnnctional relations: 

t\u)t\u + A) = 1 + t\u + 3A) 

where the off-critical transfer matrix is donble-periodic 

t\u)=t{u + Ti) ; t'{u P ine) = C ■ t'{u) ; q = (4.54) 

This means that du logt'{u) is also periodic. The prime indicates that the expressions 
correspond to the massive case whose fnnctions are determined in terms of ellip¬ 
tic thetas to differentiate them from the critical expressions used before which were 
trigonometric functions. 

Now we rewrite the functional relation as t'{u)t'{u — ^) = ^ A t'{u — , and after 

shifting u we get 

t\u + ^)t\u - ^) = 1 -h t\u) 

5 5 

Similar to the critical case, we need to kill the iVth order zeros at — ^ and poles 

at ^ by f'{u) which satishes 

+ = f{u) ■ f'{u + me) = f'{u) (4.55) 

The resulting function l'{u) will be analytical and nonzero in the required domain. So 
we dehne 


t\u) = fiufl'iu) 
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The solution for f'{u) compatible with the analytical structure is 


q) 


M'f + hP)Mf + T’P) 

MT + hp)Mf + T,p} 


(4.56) 


6 O O ' 

where the periodicity requires p = q^. Introducing the variable ^ we write 

the functional equation as 


t'{x — i—)t'{x + i—) = 1 + t'(x) (4.57) 

3 3 

The periodicity box in the variable x is 3fJe(a;) G (—vre^, ttc^) and ‘Am{x) G (—|, f), 
which is the analogue of the physical strip. 


Vacuum state 


We divide 4.57 by t and use the functional relation to write 


t'{x — i^)t'{x + i^) V{x — i^)r{x + i^) ^ ^ 


f(x) 


l'(x) 


= l + t'-Hx) 


After taking log (both sides are ANZ in the physical box) we solve it in Fourier space. 
The functions are periodic with period x = x + V, where V = ^. So we expand 
them as 


h{x) = ; 


V_ 

hk = - f" h{x)e-^^^^dx 
V J_v 


u = 


27r 

V 


and solve the equation for log/'(a;) as 


log 4 = 


log(l + T ^)fc 


our in real space 

1 

logl'(x) =—(p^iqlog(l + t'~^(x)) := — / dx(p^(x — x')log(l + t'~^(x')) (4.58) 

V . -V 
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where 


k 


^iujkx 
1 — 63 ^^^ — 


(4.59) 


Now we restore t': 


\ogt'{x) = N\ogf'{x) — log(l + t' ^(x)) 


(4.60) 


This is the ground-state massive TBA on the lattice. 

We take the continuum limit such that we scale x ^ x + logiV and N —)■ cx): 


, 37r 3i(a; + logiV) 

hm TV log/ — +---,p 

w^-oo 10 5 


In order to have a hnite limit we send p —)■ 0 and obtain 


hm N\ogf{^ + ^\ p) = 2A/3(e -1- N'^p^C) + 0{p^N^) 

v^-oo 10 5 


If we choose 


P = 


mL 1 /i 
Wl>N ^ N 


(4.61) 


This is equivalent to shifting x as x ^ x — log /r.Then after shifting the x variable we 
obtain the standard massive TBA on the lattice: 


logt/x) = mL cosh(a;) — ^ log(l-I-t' ^(a;)) (4.62) 

as in this limit e —)■ cx) and so t p. 

If we dehne 

gda:) ^ _|_ iQg jY _ iQg 

N^OO 

Then it satishes the massive periodic TBA 


e(a;) = mL cosh a; — <pa^log(l-I-e 


(4.63) 
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Excited states 

We repeat the same analysis as in the critical case and indicate the important and 
new elements in the analysis. From the numerics, we assume the existence of single 
real zeros and complex conjugated short strings. However, contrary to the critical 
case we do not assume that they are located exactly at ±^. The relevant functions 
to eliminate these zeros with the required periodicity are 


a()(a;,r) 

a[{x,p) 


M^,r) 





(4.64) 


(Jo(a;,r) has a single zero at a; = 0 within the periodicity box, while a[{x,q) has one 
at and another at —f|. So eliminating singularities on the upper half plane at /3i 
> 0) we use cr[(x — A — y,p) while on the lower half plane at /dj (!fle(/3i) < 0) 
we use a[(x — /3i -h p). A short two string with locations f3k = Uk L ivk can be 
eliminated by 


a[(x -Uk - — - ivk,p)(T[{x -Uk + ivk + —,p) = o'^ix - Uk,Vy,p) (4.65) 
6 6 

For a state with one-strings at x^ and short two-strings at xf = {uf,vf), the calcu¬ 
lation analogous to the massless case leads to 


logF(a;) = N\ogf{x) + J2± logK(a^ - > ^)) + Ey± ^ogia'^ix - uf, vf, p)) 

-ip,^\og{ll^a'Q{x - x^,r){l + t'-^{x))) 


(4.66) 


We do the continuum scaling as we did in the vacuum case while, in addition, we 
suppose that the one string and the two strings scale as Xk ^ Xk T log “ , then we 
have to take into account that in the scaling limit 
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lim ao(a; ± log f - x^, r) ^ = ao{x - aj) 
lim a'^ix ± log ^ - Uk,v^,p) = 0-2(0; - Uk,v^) 

N^oo ^ 

where we retain the trigonometric expressions (Tq{x) = tanh(^) and 0-1(0;) 
-^) tan(y + x^), and dehne 


(4.67) 
tan(y — 


0 - 2 ( 0 ; - Uk,v^) = 0 - 1 ( 0 ; -Uk - — - ivk)cri{x -Uk + ivk + —) 

6 6 

we can dehne the function 

^40 ^ _ XQ)~^t'{x + log —) 

N^oo jj, yU 

and it satishes the equation 

e(o;) = mL cosho; + log 0 - 2 ( 0 ; — Uj, Vj) — ^p-k log(o-o(o; — aj) + (4.68) 

j 

The parameters as before, satisfy the following equations = —1. 

This is the massive excited state for the Lee-Yang model with periodic boundary 
conditions. 

4.3.2.2 Seam 

In a similar derivation to the critical periodic model with a seam, we derive the massive 
one here. The periodic transfer matrix in presence of a seam satishes the following 
relation; 


TAu)TAu + X) 


s{u-\-X + 0sU-u-0 f s{u-\- X)s{X- u)\^ \ s(-«-hC) f s{u)\^ ^ i 

-- \ — 


with 
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, _ s{X)s{u + 2A + ,^) / s(A)g(-u + 2A) ^t( 

s(m + A + .^)s(m + 3A + .^) \s(m + A)s(m + 3A)/ ^ 

we obtain the relation 


t'{u)t'(u + A) = 1 + t' (u + 3A) 


Then, writing t'{u) = f'^ ^{u)f'{u + and killing the zeroes of the normaliza¬ 
tion, exactly as we did before and making the change of variables n = -1- we 


obtain 


fix) = 


Mf + ^,p)Mf + ^,p) 


(4.69) 


and the periodicity requires that p = qs as before. 


Vacuum state 


Since this state has no zeros in the physical strip we get that: 


logf'(a;) = (iV-l)log/'(a;)+ log/'(a;--(y9,^log(l + f' ^(x)) (4.70) 


where the convolution * and were defined in (4.58) and (4.59) respectively. This 


is the massive ground-state TBA on the lattice with a seam. 

We take now the scaling limit. We calculated the scalings of those functions before in 
the previous section. Only log/'(x — ^) was not considered. As x —>• x -|- log ^ and 
taking ^ = —i|(^ -f log ^), we get that the continuum limit of /'(x — ^) is exactly 
its trigonometric form /(x — ^) where 


ix dvr ix Tn 
/(x) = - tan(y + —) tan(— + —, 


Defining = lim t'(x -|- log —) we get 
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e{x) = mL cosh X + log f {x — — if-k \og{le (4.71) 

In other words, the ground-state massive seam TBA becomes: 


e(a;) = mL coshx-I-logT_(-^-I-x) — 9 ? ^ log(l-I-e (4.72) 

Excited states: 

To eliminate the zeros of the 1-string and short 2-strings with the required periodicity, 
we use ctq^x, r) and a'i{x,p). As before a short two string with locations (3k = UkEivk 
can be eliminated by 

aRx-Uk - — -iVxkMDA^-^kP^v^kP yr,?) = o^KX - Uk,v^,q) (4.73) 
0 0 


where vi and V 2 are not necessarily equal. 

Now suppose in the state we have one string at Xq (either in the upper or in the lower 
half plane) and short two strings at x^ = A calculation analogous to the 

massless case leads to 


logT(a;) = (iV- l)log/'(a;) + log/'(a;-0 + ^^log(ao(a;-a;^,r)) + 
^og{(T'^{x - u^,vf,p)) -ip-k log(n± (^o{x - x^, r)(l + t'-^a:))) 


(4.74) 


We already know the scaling limit of this relation. Dehning the function 


ehU = lim (To(a; ± log — — Xn ) ^t'(a;-|-log —) 

N^oo p, p 


(4.75) 


It satishes the equation 


e{x) = mL cosh a; -|- log/r(a; — f) + log a 2 {x — uj, vj) — ip-k log((Jo(a; — Oij) P e 
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or using the relation derived before /(x — .^) = T_(y + x) ; b = 3 + ^ then 


€{x) = mL cosh a; + logT_(— + x) + \oga 2 {x — Uj, Vj) — 99T*clog((Jo(a; — a) + e 

(4.76) 

which is the massive excited-state TBA of the Lee-Yang model with a seam. 


4.3.2.3 Boundary 

Similar to the massless case we dehne: 


t{u) = Ss{u) 


{}1{2u — A, q) 


ili{2u -|- A, g)-di(2n — 3A, q) L-di(n -|- A, q)'&i{u + 3A, q)i 


h'i(A,g)h'i(n -h 2\,q) 


N 


D{u) 


with 


= 


1 , 


_ '&l{\,qY'dl{u+^+2\,q)'&i{u—i+2\,q) _ 

'i9i (ti+^+A,ij)i?i (ii+5+3A, 5)191 (ti-^+3A,5)i9i (ti-^+A,^) 


S = 1 


S = 2 


Then the normalized transfer matrix satishes the universal TBA functional equation 


t{u)t{u -|- A) = 1 -|- t{u + 3A) 


We write 

t'{u) = nuf^g'{u)l'{u) 


To kill the Yth order zeroes and poles, and after shifting n, and with p = q^we have 


^i(f+ ff, p)Mf + ^, p) 
^2(f+ ff, p)Mf + ^, P) 
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(r,s) = (l,l) sector: 

We also want to eliminate the zeroes of order 1 which originate from the normalization 
and from the analytical strip. In a similar analysis of the critical case, and taking 
into account the periodicity implied by the normalization, we deduce that the factor 
that we need to insert and compatible with the analytical structure is: 


9'i,L 


u = 


- l,pi) 

i(f - - |,p3K(f - f.pi) 


In terms of the x variable we get: 


9'i,L 


X) = 


,p3)i?,(f+i,pj)pi(f - i.p5)tf2(f - +s.p") 


+ l,p-^)U'-i - i.p=)A(f - ■S.P=)!’i(f + r^.P'^) 


The priodicity of g is half that of /, due to the fact that the order 1 normalization 
is a function of 2u, rather than u, hence the period is divided by two, or equivalently 
the nome oi g is p^. 

g[ satishes the relation 


9[,i{x - i^)g[,i{x + i^) = g[^i{x) (4.77) 

We eliminate the zeros of the 1-string and short 2-strings with the required periodicity, 
we use (Tg(a;,r) and a 2 {x — Uk,Vv,p). As before a short two string with locations 
/3k = Uk E ivk can be eliminated by 

ZTT ZTl 

a[{x-Uk - — -ivik,q)(T[{x-Uk + iv 2 k + yr,?) = a^ix - Uk,v^,q) (4.78) 

D D 

where Vi and V 2 are not necessarily equal. 

In the X = §i{u — ||) variable they are located at Xq = ^ = logA^ -|- Q;=*"and at 
x^ = -|- log 3- /3f. As we would like to take logarithm we need functions free of 
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zeros and poles on the physical strip. Those functions were dehned before as aQ{x,r) 
and a 2 {x,p). 

With these functions we parametrize the normalized transfer matrix eigenvalue as 

M 

= f{xf^gl^{x) JJ a^ix - r) JJ - Uk, v^, q)l'{x) (4.79) 

± j=i 

which ensures that l'{x) is ANZ in the physical strip. 

We hnd: 


log/'(a;) = -(pe*log(JJcro(x-a:o)(l + ^' ^(a^))) 


where the convolution -k was dehned in equation (4.58). Restoring t we obtain 


logf'(x) = 2N\ogf'{x) + log^^(x) + log(cT(,(x - Xq )) + J2j,± ^og{a' 2 {x 

-Te * log(n± - 3^0 )(1 + 

The parameters are determined by t{x)\^^^.^in = —1. 

In the scaling limit, g[^i{x) —)■ 5 'i,i(x) = R7^(y + x)i?7^(y — x). 

Using 

^ ^\Qg^ — XQ)~^t'{x + log —) 

N^oo g, g 


Uk,Vy,q)) 


(4.80) 


(4.81) 


The parameters satisfy e'^^®Vo(x — = — 1. 

And e(x) satishes the equation 


e(x) = 2mL coshx — logi?i(y + x)i?i(y — x) + Xli log'^ 2 (^ — Uk,Vy,q) 

^ (4.82) 

— if -k log((To(x — «='')+ 

This is the general massive TEA in the (1,1) sector. It is important to mention that 
for the groundstate of this sector, there are no short 2-strings, hence ai doesn’t appear 
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and the ground state massive TBA of this sector on the lattice is: 


\ogt'{x) = 2N\ogf{x) + \ogg[^^{x) + J2^\og{aQ{x-x^)) 
-if, ^ log(n± (^o{x - )(1 + 

while the scaled ground-state massive TBA is 


(4.83) 


e(a;) = 2mLcosha; —logi?i(-^— x)-\ — ip-k\og{(7Q{x — Ci^) + e (4.84) 


(r,s) = (l,2) sector 

Similar analysis follows in this sector, where we have to take into account the contri¬ 
bution of Ssiu). 

For s = 2 , 5*2 (m) is proportional to f {u + ^) f {u — , therefore we should account 
for their respective zeroes and poles. 

The other contributions from the normalization and the analytic strip are the same 
as before, and they were included in g[ i(m). Then we need to dehne the new g[ ^ term 
as: 


X f{u + Of'{u - 0 

In terms of the x variable we get: 


^' 1 , 2 ( 2 ^) 




ix 

2 

ix 

2 


,phMf+fxhMf-^,phM 

xfix-f)fix+f) 


ix 

2 

ix 

2 


r2xhMf 






(4.85) 


digix) satishes the relation 


9'l,2i.X 




9'i,2{x) 


The zeros short 2-strings are eliminated by ^^(a; — Uk,Vy,p) 

In the X = §-i{u — ^) variable they are located at = ±f| -|- logiV + As we 
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would like to take logarithm we need functions free of zeros and poles on the physical 
strip. We take a[{x,p) to accommodate for the zeroes of the short strings. With these 
functions we parametrize the normalized transfer matrix eigenvalue as 

M 

- Uk,V^,p)l\x) 
i=i 

then we obtain 

logt'(a;) = 2N log f'{x) + logg'^ ^ix) + ^ log(a 2 (a; - Uk, v^,p)) - log(l + t'~^{x)) 

(4.86) 

The parameters are determined by t '= —1. 

In the scaling (iV —)■ cx)) limit, as a; —)■ a; + log ^ and ^ —)■ —i^{^ + log then 

ITT ?7r 

g[^,{x) ^ g,,,{x) = R-/{- + x)R^\- - x) 

for 6 = 3 + ^ in the expression of R^. 

Using = limAr_,.oo t{x + log ^), it satishes the equation 


e(a;) = 2mLcosha; - logi?. + x)R^ Y ~ ^) + Zli log<^ 2 ( 2 ^ ~ Uk,Vv,p) 

—(p ic log(l + 

This is the general massive TEA in the (1,2) sector. For the groundstate of this sector, 
there are no short 2-strings, hence ai doesn’t appear and the ground state massless 
TEA of this sector on the lattice is; 


logt'(a;) = 2iVlog/'(a;) + log^^ 2 (a^) + -<^*log(l + t' ^(a;)) 
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while the scaled ground-state massive boundary TBA in the (1,2) sector is 

e{x) = 2mLcosha; — log-R^^(^ -1- x)R^^{^ — x) -\—ip ^ log(l -f- (4.89) 

We have derived all the TBAs for the periodical boundary model with and without 
a seam, and for the boundary model in both critical and massive regimes for ground 
and excited states. 



Chapter 5 


Energy corrections: from Lee-Yang to 
SYM 


As mentioned in the introduction, AdS/CFT correspondence jT] relates the string en¬ 
ergies on the AdS^ x background to the anomalous dimensions of gauge invariant 
operators in maximally supersymmetric four dimensional SU{N) gauge theory. The 
correspondence is particularly useful in the large N limit when it can be described 
by a two dimensional integrable held theory |9H]. The advantage of the integrable 
two dimensional point of view lies in its non-perturbative nature, which nevertheless, 
can be used to calculate perturbative quantities, too. Indeed, the leading hnite size 
correction of a two particle state can be described in terms of the asymptotic Bethe 
ansatz and the generalized Liischer formulas |99], which, when expanded in the cou¬ 
pling constant, provides the exact perturbative anomalous dimension of the Konishi 
operator up to seven loops. 

To show this resemblance, we start by calculating higher order corrections for the 
energy of a one particle state in the Lee-Yang model. We do this because it turns 
out that the F — term corrections can be calculated in an analogous way for several 
other models. We deduce a general result that can be applied to integrable models 
including the AdS superstring sigma model. 


94 
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In the bulk Lee-Yang model, we know that the multiparticle state energy is given by 


N 

E^nj}{L) = cosh^j — m 
j=i 

The ground state TEA is given by 


+00 


— cosh6'log(l-I-e 


e{6) = mL cosh^ 



6^') log(l-1-e 


(5.1) 


where (j){6) = —idg log S{6) and S{6) is the scattering matrix of the model |6ni| 
In the large volume limit, it reduces to: 


e{d) = mL cosh6 

Substituting this into the energy equation we get that the leading hnite size correction 
is given by: 


/ -I-OO jg 

— cosh01og(l + e-™^“"^®) 

ZTT 


^+00 JO 

= -m —cosh 


' —CXD 


271 


(5.2) 


while for a moving one particle state TEA we need to do one more iteration. The one 
particle TEA equation is 


e{9) = mL cosh 9 + log 


5(0 - 9o) 
S{9 + 9o) 


1 



— 00 


O') 


S{9 -|- 9o) raL cosh S', 

5(0 - 0o) ^ 


(5.3) 


The 00 is not arbitrary. It is determined from the relation e(0o) = i7r(2n-f 1) therefore 
(neglecting the third term for large L) we have 
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mL cosh 9o + in — log S'(26'o) = in{2n + 1) 


or 


mL cosh. 6 q — logS'(26'o) = 2inn 

For a standing particle, n = 0 thus the hrst two terms need to cancel each other. For 
large L this can only happen if 9q causes a singularity in S{29q), hence for 0o = y- 
Continuing the analysis for a standing particle, the contribution to the F — terms 
will come from the energy correction 

+00 

E,{L) = -m [ dgcoshg ^l^ + y e— 

—oo 

Hence the hrst order correction for a standing particle state if given by: 

+ 00 

E,{L) - Eo{L) = -m J d9 cosh9{^^^^ - i)e-'^Lcoshe 

— OO 

Now using that 0o = ^ + id, where 5 is a small correction , and the fusion relation of 
the Lee-Yang S{9+ in/3)S{9 — in/3) = S{9) we get the energy correction for the one- 
particle state. 6 will lead to /i — term corrections whereas the F — term contribution 
will be: 


+ 00 

Ei(L) - Eo{L) = -m J d0cosh0(^(y -9)- 

— OO 

This hnal result is general. It applies to the F — term contribution of the Lee-Yang 
[60| . the sinh-Gordon |99| . some sectors of the sine-Gordon and several other models. 
It is conjectured to be the relation that gives the higher order F correction terms for 
integrable model TBAs. It is generalized to the AdS model. 

With the AdS/CFT conjecture and the identical 4-loop corrections that we hnd on 
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both sides of the duality, we can calculate higher order corrections for the superstring 
model and assume that those results will hold on the other side of the duality.The 
direct perturbative gauge theory calculations are very cumbersome and have been 
evaluated for the Konishi operator up to 4 CL^d 5 loops only [Ton [TOO] , where the 
wrapping part of the correction matches exactly the Liischer type hnite size correction 
|99j . Thus the integrability based “string theory” techniques provide a way to go 
beyond the available perturbative calculations and collect information about higher 
order results. 

The /9-deformed theory is an = 1 supersymmetric, exactly marginal deformation of 
the A/" = 4 SYM and has been always the testing ground of the AdS/CFT duality. The 
dispersion relation of the excitations was calculated exactly in [104] . which together 
with wrapping corrections [1051 1106] provides the exact anomalous dimension of a 
single impurity operator Tr(XZ^). On the string theory side this operator corresponds 
to a one particle state in the su{2) sector in hnite volume. The leading hnite size ehect 
follows from the momentum quantization, or asymptotic Bethe Ansatz equation [Ml. 
As the volume decreases Liischer/wrapping type correction becomes important [991 
110211108] and for an exact description they have to be summed up. An educated way 
to take all hnite size correction into account would be to extend the Thermodynamic 
Bethe Ansatz of the A/" = 4 theory [1091 111111112111131111411115] to the /3-deformed 

case. This has been recently done in ina. 

We are interested in the 13 -deformed theory at a particular value of f3, namely for 
the simplest nontrivial 13 = \ and evaluate the Liischer correction at NNLO in the 
coupling. This value of f3 is interesting as it correspond to a particular orbifold theory, 
too. As it was conjectured in [118] the anomalous dimension of the su{2) particle 
(Tr(XZ)) in the (3 = 1-deformed theory coincides with the anomalous dimension of 
the si(2) particle {Tt(DZZ)) of the orbifold theory. 

In the results, multiple zeta functions appear. Multiple zeta values (MZVs) are the 
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generalization of zeta fnnctions; 


C(®1) ®25 • • • 5 ®n) ^ ^ 


jl>i2>->in>0 


Jl J2 * * * Jn 


and are trivially related to the values of nested harmonic sums taken at inhnity. A 
recent review is presented in |123) . 

We calculate the anomalous dimension of the operator Tt(XZ) in the P = \ deformed 
theory. It is a very special state as it has a vanishing rapidity n = 0, (p = tt), which 
is not effected by hnite size corrections. In contrast, the energy is shifted by vacuum 
polarization effects as 

E = Eaba + Ai? (b-5) 

The asymptotic Bethe Ansatz energy Eaba is simply the dispersion relation of a 
standing particle: 


Eaba = = a/1 + 16g‘^ 

= l + 8g^ - 32g^ + 256p® - 2560/ + 28672p^° - 34406%^^ ^ q 


while AE corresponds to the wrapping interactions and has the expansion 


AE = AE^g^ + AE^g^^ + AEeg^'^ + ... 


,12 


(5,6) 


In |in2) the LO correction was calculated 


AE, = 128(4C(3) - 5C(5)) 


(5,7) 


The NLO correction turned out to be im : 


ABs = -128(12C(3)^ + 32C(3) + 40C(5) - 105C(7)) 


(5.8) 
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And in |128] we evaluated the NNLO correction and obtained 


AEe = -128(48C(3)2 - 592C(5) - 24C(3)(8 + 15C(5)) - 322C(7) + 1701C(9)) (5.9) 


which completes the anomalous dimension of the operator up to six loops. 



Chapter 6 


Conclusions 


In the first part of this thesis I presented an introductory review of the Lee Yang model. 
The in chapter 3,1 initiated the form factor program for purely transmitting integrable 
defect theories. I restricted for a single particle type, but the extension of the program 
for diagonal bulk scatterings and diagonal transmissions is straightforward |92j . We 
laid down axioms for the form factors of operators localized both in the bulk and also 
on the defect. We determined the solutions of the consistency requirements for bulk 
operators in terms of the bulk form factors together with the transmission matrix. 
These form factors determine the correlation functions of bulk operators, which we 
elaborated in details for the two point functions. In the case of defect operators we 
gave the general form of the solutions and explicitly calculated for the free boson and 
for some operator in the Lee-Yang model. We also described how the fusion method 
can be used to generate new form factor solutions from old ones. 

In the analysis of the Lee-Yang model we observed relations between the defect oper¬ 
ators which should have the origin in defect conditions. The lack of the Lagrangian 
dehnition of the model prevented us to analyze this question. In order to achieve this 
aim one has to analyze the simultaneous integrable defect and bulk perturbations 
of the defect Lee-Yang model using conformal perturbation theory and establish the 
relation between the bulk and defect couplings, which maintains integrability. This 
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approach then can be used to derive defect conditions which will provide relations 
between helds living on the defect. In pushing forward this program one has to solve 
the defect Lee-Yang model hrst. The explicit knowledge of the correlation functions 
together with the structure constants will make it possible to bridge the operators ap¬ 
pearing in the form factor program to their UV counterparts by analyzing the short 
distance behavior of the two point functions obtained from the explicit form factor 
solutions. 

Another interesting problem is to see that we have as many polynomial solution of 
the form factor equation as many local operators existing in the theory. We have 
seen this coincidence in the case of the free boson. In the case of the Lee-Yang model 
the nontrivial mixing between the left and right degrees of freedom and the various 
cancellations between the leading order scaling terms prevented us from performing 
this analysis. Possibly a more careful analysis along the line of ra would clear up 
this point as well. 

We have analyzed the free boson and the Lee-Yang model, the method, however, 
has a straightforward application for the sinh-Gordon model adopting ideas from the 
boundary form factor solutions [93l IM] . 

The defect form factors in the Lee-Yang model can be tested by extending them for 
hnite volume and comparing to direct TCSA data. They also can be used to build 
up hnite temperature defect correlation functions. These are direct generalizations of 
the related boundary analysis developed in j95l 1^ . 

Whereas in the second part, in chapter 4, we dehned the Lee Yang model on the lattice 
as a special case of the Al Forrester-Baxter model. We dehned the face weights, for 
both massless and massive models and then dehned the transfer matrices for the 
periodic boundary conditions, with and without a seam and also the double row 
transfer matrices for the boundary case. Analyzing their analytic structure by hnding 
the zero eigenvalues of those matrices using numerical methods we were able to provide 
a systematic classihcation and deduce their hnitized characters. We were also able 
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to deduce a direct bijection between the paths, the zero structure and the Virasoro 
modes. Given the analytic structure of the transfer matrix for any energy eigenvalue, 
we can associate it in a one-to-one correspondence with a state in the Virasoro module. 
For the boundary case, I classihed the boundary flows between the states of the two 
modules of the model exactly. More interestingly, we applied this to the case of a seam 
where the operator content and the flow patterns are much richer, and the method 
accurately describes all the flows in terms of the states, the zero structure and the 
path diagrams. 

With the knowledge of the analytic structure, and using the functional relations of the 
model in the three different boundary conditions, one can solve the TEA equations for 
all cases in the massless and the massive regimes. The seam TEA is also consistent 
with the transmission matrix suggested in the bootstrap. The results also conhrm 
what was found before through the TEA approach. However, before, the analytic 
structure was guessed. On the lattice, in contrast, it is exactly determined. 

The methods of this part are general in the sense that they can be applied for a variety 
of models of more complexity. A future line of research can work on solving other 
integrable models. 

In the large volume limit, I calculated the Luscher correction terms of the TEA 
energy of a one-particle state to higher orders. This generalization is important as 
similar techniques are used in the analysis of energy corrections due to wrappings for 
superstrings in the sigma model. I display those ideas in the end of the dissertation 
as this is an important direction of work in string theory models, particularly in the 
AdS/CFT correspondence. 



Appendix 


Calculation of the anomalous dimension from AdS/CFT 

In this appendix we explain how to obtain the leading Liischer correction for a one 
particle state in the (^ = \ deformed theory by recalling the available formulas from 
the literature. 

As the anomalous dimension of the single impurity operator corresponds to a standing 
particle state, we have to calculate the standard Liischer correction. It describes how 
the one particle energy (dispersion relation) 

E = E{p = tt) = + 16g^ sin^(^) = a/1 + 16g^ = e (6.1) 

is modihed due to vacuum polarization effects. The p = tt momentum corresponds to 
vanishing rapidity u = 0 which is protected by symmetry. Consequently the leading 
hnite size correction is the energy correction which originates from virtual particles 
propagating around the circle and can be written as 

AB = - f; f ^sTr(S«;(,,0))e-‘«<'»^ + 0(!,“) (6.2) 

Here we sum up for all bound-states of charge Q of the mirror model, whose momenta 
are q, Sq^ describes how they scatter on the fundamental particle and e(q) denotes 
their mirror energy. We have to expand this expression in NNLO in g. The above 
form exactly describes the energy correction up to the order g^^ only since at the order 
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double wrapping effects will contribute, too. Let us analyze the g dependence of 
the various terms. 


Exponential factor 

The mirror energy has the following parametrization 


-^(,) ^ ^ {q,Q) 




q + iQ 

^9 



(6.3) 


Matrix part of the scattering matrix 

The contribution of the scattering matrix can be factored as 


sTr(^g(g, 0)) = S^caiariq, 0)sTr(5^“i2.^(g, 0))^ (6.4) 


where the super-trace of the matrix part contain the contributions of all polarizations. 
The various polarizations of the mirror bound-states can be labeled in the super-space 
formalism as 9^) see 1124) for the 

details. The super-trace in the /3- deformed theory evaluates in the previous basis as 

Q Q-2 Q-l Q-l 

sTr(^ra?ix(^> a;)) = ^ SBlj{z, a;)+^ SB2j{z, x)+i ^ SFlj{z, q)-i ^ SF2j{z, q) 

j=0 j=0 j=0 j=0 

(6,5) 

The S-matrix elements can be extracted from |99| . In calculating the corrections for 
the su{2) representative neither of the S-matrix contributions depends on j and they 
read explicitly as 


SBlj{z, x) 


z^ — fji 

z~ — x+ gi 


SB2j{z, x) 


Z~^ — X (1 — X^Z ) X fji 
z~ — a;+ (1 — x~z~) a;+ gi 



2 


Z~^ — X 1]i ?72 
— a;+ gi r/2 


z~^ — (1 — x~^z ) X ?7i r/2 

(1 — x~z~) a;+ r /2 


SFlj[z, x) 


SF2j{z,x) 
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The appearing string frame factors can be written as ^ = y ^ and • The 

parameters depend on the momentum the usual way 


^ cotfii 
= —Jg — 


+ 16^2 51^2 


( 6 . 6 ) 


which in our case results in 


I 1 , s X 

= ±x = ± —(1 + e) ; ^ =-1 

Ag^ x+ 


With these variables the super-trace of the matrix part takes a particularly simple 
form 

0-T(^ 4- 9/0 (T — 

(6.7) 


sTrfc-( 2 ) . .. ^ 2x{l + 2Qz {x - z+)) 

^^n^matrix[UX)) (j; - z-){l + XZ-) 


We have checked that the matrix part of the sl{2) representative of the orbifold 
model gives the same result. There, the deformation is such that the fermions do not 
contribute as the two undeformed fermionic S matrix elements are the same |1251lin2] . 


Scalar part of the scattering matrix 


The scalar part of the scattering matrix of a charge Q bound-state can be obtained by 
multiplying the scalar factors of its individual scattering constituents. The charge Q 
bound-state composed of elementary magnons as z = {zi,..., zq), such that 2 ;“ = Zi 
and Zq = z+and the bound-state condition is also satisfied z^' = zyy^. Thus the full 
scalar factor as the product of the elementary scalar factors turns out to be |99j : 


Q 


Sscalar y x'j * 5 * 


sl{2). 


1=1 


ZiyX) = lie 
i=l 


Q 

n 


— 2ia{zi,x) 


X 


zj — x~ 1 — 


( 6 , 8 ) 


In calculating the Liischer correction we have to evaluate this expression when 2 ; is 
in the mirror kinematics (|^“| < 1, \z~^\ > 1). The analytical continuation has been 
carefully elaborated in |126j : 
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SscLlari^^^) = (6.9) 

A particularly good feature of the formula is that both expressions depend on 2 ;^ only 
and not on the individual zf. Explicitly the su{2) scalar factors read as: 

(2+-I-)(z+-X+)(l-^)(l-^) 

while for the case |x^| > 1 following HSl we can write 


—ilogSQ^i(2;,a;) = ^{z'^,x'^) — ^{z'^,x ) — $( 2 ; + ,x ) 

+ ^ [—'I/{z~^,x~^) + 'I/{z~^,x~) — 'I/{z~,x^) + \l'(^“,a;“)] 

1 r {z^ - x^){x-- - 

We have to be careful as the conventions of ra are different from ours. To turn 
into our conventions one has to replace 2 ;^ —>■ z^ and x^ —>■ There are integral 

representations for <h andT as 


*h(xi,X2) = i 


dwi f dw2 


27ri 2m wi — X 1 W 2 — X 2 


log 


T{l + ig{wi+w^ -W 2 -W 2 )) 
r(l - ig{wi + -W 2 - W 2 ^)) 

( 6 . 11 ) 


T(xi,X2) = i 


dwo 


log 


r(l + i 5 ((xi + - m;2 -'W^2 )) 


( 6 . 12 ) 


Tc^ 2m W 2 - X 2 ° r(l - ig{xi + x^^ - W 2 - w^^)) 
where the integrations are for the unit circle. They are well-dehned provided none of 
the X, lies on the unit circle. 


As by now we have collected all the necessary formulas we turn to analyze their weak 
coupling expansions. 
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Weak coupling expansion 


Our aim is to calculate the weak coupling expansion of AE for L = 1. In doing so we 
decompose the integrand of the Liischer correction 


/ ^P(q,Q)Uq.Q) 


into a simpler rational part 


4x^(1+ 2Qz (x-z+)f /z y 

((xz-)^ - l)(x^ - (z+)^) J 


(6.13) 


(6.14) 


which contains both the matrix part and the rational part of the scalar factor, and 
into the more complicated S part: 


ilogE(q,Q) = 2(4>(z~^,x~^) — $(z~^,x ) — $(z ,x~^) + $(z ,x )) 

—'l/(; 2 ^, a;’*') + \['( 2 :’'', a;“) — 'L(z~,x~^) + 'L(z~,x~) 

We expand these functions in as 

P(g, Q) = Ps{q, Q)g^ + Pio(g, Q)g^^ + Pi 2 (g, Q)g^^ + ■■■ (6.15) 

S(g, g) = 1 + S 2 (g, g)/ + S4(g, Q)g^ + ... 

The expansion of the rational part is quite straightforward and we obtain 

4096g^(-l + g^ + g^)^ 

(g2 + g2)4(g4 + (_1 + g2)2 + 2q^{l + g^)) 


Pioig-iQ) _ 8(7g^ + 3(g^ — 1)^ + i0g^(i + g^)) 

Psiq^Q) ~ “(^^+gW+Fi+W+W+W 


(6.16) 
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Pi 2 {q,Q) _ 2(27 + 268g2+ 704g4) + 2q^Q) 

Pw{q,Q) ~ (3 + 28g2+6V)(<j2+Q2) + (l+4g2)(l+^2_2Q + Q2) 

16(3g2 + 2q^Q) 8(-81g2 - lHq'^ + IQq^ + 27q'^Q^ + 16q^Q^) 

“(1 +V)(l + q2 + 2Q + Q2) “ (3 + 16g2)(3 + 10*72 + 7g4 _ + l0q2Q2 + 3Q4) 


In expanding the \I/ and $ functions we use the same method we used in (7). The 
expansion of the 'I/{xi,X 2 ) functions for \x 2 \ > 1 (string region) reads as follows 


'^{xi,X2) = - — {'I7{l-ig{xi + x^^)) + 'lf{l + ig{xi+x^^))) 

X2 

-^(^ 1(1 - ig{xi + - Ti(l + ig{xi + 

g3 

+ 7 ^('J' 2(1 - ig{xi + x^ ^)) + ^ 2(1 + ig{xi + ^))) + ... 

^X2 

where = (^)"'(log(r(a;)) are the standard polygamma functions. If |a;i| > 

1 then $(a;i,a; 2 ) starts at g^. In the opposite case using the identity <h(a;i,a; 2 ) = 
$( 0 , 0 : 2 ) — $(a: 7 ^,a: 2 ), being valid if |a:i| 7 ^ 1 , we can calculate the leading expansion 
of $ as 

$(0,a:) = -(7 e^-3C(3)/ + ...) 

X 

Using functional identities valid for integer Q we obtained 


Mq.Q) 


TTqs “ ® 0 1 —2 —) I—2— 


and 


^4{q,Q) — ^^2iq,Q)^ + 


64Q(l 




(g2 + g2 


32iq ( ^ fQ - iq-2 

q^ + Q- 


H 


32 


U( 


Q — iq — 2 


)-s,{ 


Q iq — 2 

2 J 2 

Q - iq-2\ fQ + iq-2 


+ 18C(3) 
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where Sn{x) are the analytical continuation of the harmonic sums Sn{N) = 
They are related to the polygamma functions a^ 

T„(Q) = (-l)"+'n!(C(n + 1) - Sn+i{Q - 1) 

Further details are given in |128). 

^For n = 0 one has to replace C(l) with je- 
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